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ABSTRACT 


Ocean  internal  waves  induce  magnetic  fields  by  virtue  of 
motion  of  conducting  sea  water  relative  to  the  geomagnetic  field. 
Measurements  of  such  Internal  wave  Induced  magnetic  fields  by 
sensitive  magnetic  instruments  below  the  ocean  surface  require 
use  of  a  protective  enclosure  or  buoy.  In  this  paper  an  anal¬ 
ysis  is  presented  of  the  effects  of  the  enclosure  on  the  mag¬ 
netic  fields  and  their  spatial  gradients  when  the  enclosure 
(buoy)  is  spherical. 

Explicit  analytical  results  are  presented  for  the  magnetic 
field  and  the  magnetic  field  gradients  for  the  case  in  which 
the  sphere  radius  is  small  both  in  relation  to  the  spatial 
period  of  the  Internal  wave  field  and  in  relation  to  the  depth 
at  which  the  buoy  is  submerged.  In  addition,  analytical  results 
are  presented  for  the  temporal  spectra  of  magnetic  field  gradi¬ 
ents  as  these  would  be  measured  within  the  buoy  when  towed  at 
speeds  substantially  greater  than  the  maximum  group  velocity  of 
the  Internal  wave  field. 
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I.  INTRODUCTION  AND  SUMMARY 


The  measurement  of  magnetic  fields  and  magnetic  field  gra¬ 
dients  In  the  ocean  that  are  Induced  by  the  relative  motion  of 
sea  water  with  respect  to  the  geomagnetic  field  typically  re¬ 
quires  that  the  measuring  Instrument  be  enclosed  In  a  measure¬ 
ment  chamber  or  buoy.  For  example,  when  the  Instrument  Is  a 
superconducting  gradlometer  (SQUID  sensor)  the  minimal  enclosure 
would  be  the  cryogenic  unit  Itself.  The  major  contributors  to 
the  Induced  magnetic  fields  are  believed  to  be  ocean  surface 
waves  and  Internal  waves,  the  latter  providing  the  dominant  con¬ 
tribution  at  frequencies  below  about  1  mHz  for  a  moored  sensor. 
Vfhereas  formulas  for  computing  the  Induced  magnetic  fields  In 
the  ocean  are  available  [1],  no  analysis  appears  to  have  been 
carried  out  for  estimating  the  effects  of  the  enclosure  on  the 
detected  magnetic  field. 

Basically,  the  buoy  would  Introduce  two  types  of  disturb¬ 
ances:  fluctuations  of  the  magnetic  field  arising  from  the 
localized  flow  pattern  dictated  by  the  hydrodynamic  properties 
of  the  buoy,  and  perturbations  due  to  the  sea  water-air  discon¬ 
tinuity  at  the  boundary  of  the  enclosure.  The  latter  effect 
arises  from  the  requirement  that  the  normal  component  of  electric 
current  as  well  as  the  normal  component  of  velocity  vanish  at 
the  buoy  surface. 

The  Induced  magnetic  field  within  the  buoy  Is  determined 
essentially  by  an  Integral  taken  over  the  entire  volume  occupied 
by  the  velocity  field  exterior  to  the  buoy.  If,  as  would  nor¬ 
mally  be  the  case,  the  volume  occupied  by  the  Internal  wave 
field  Is  much  larger  than  that  occupied  by  the  local  flow  around 


the  buoy,  it  is  reasonable  to  suppose  that  the  contribution  to 
the  induced  magnetic  field  from  the  local  flow  velocity  field 
may  be  neglected. 

A  more  severe  perturbation  would  be  that  due  to  the  discon¬ 
tinuity  at  the  buoy  bounding  surface.  One  readily  convinces 
himself  of  the  Importance  of  this  effect  if  one  recalls  that 
the  asymmetry  of  the  gradients  of  the  magnetic  field  in  the 
water  is  removed  once  the  measurement  is  carried  out  in  air 
(i.e.,  within  the  buoy):  rotating  the  axis  of  a  hypothetical 
SQUID  sensor  in  air  by  90  deg  would  produce  no  change  in  the 
measured  gradient,  whereas  different  results  would  be  obtained 
if  such  a  rotation  were  carried  out  in  sea  water.  This  is 
simply  a  consequence  of  the  fact  that  the  curl  of  the  Induced 
magnetic  field  vanishes  identically  in  air,  but  not  in  sea 
water,  where  there  exist  localized  electric  currents. 

In  this  paper  an  analysis  is  presented  of  the  effects  of 
the  discontinuity  at  the  buoy  walls  on  the  magnetic  fields  and 
gradients  Induced  by  Internal  waves  in  the  ocean.  An  exact 
solution  for  the  fields  within  a  buoy  of  arbitrary  shape  appears 
rather  difficult.  Here,  for  simplicity,  the  buoy  is  modeled  as 
a  sphere.  Although  this  shape  hardly  constitutes  a  practical 
contour  for  a  submerged  Instrument  package,  it  does  provide  a 
fairly  tractable  model  for  estimating  the  general  trend  of  the 
boundary  effects.  In  addition,  the  methodology  can  provide 
guidance  for  future  analyses  of  more  realistic  shapes.  Even 
in  the  case  of  the  sphere,  an  exact  solution  is  not  trivial. 

The  approximations  underlying  the  present  analysis  are  that  the 
sphere  radius  is  small  by  comparison  with  the  spatial  wavelength 
of  the  unperturbed  Internal  wave  field  and  by  comparison  with 
the  depth  of  the  buoy  below  the  ocean  surface. 

The  first  of  these  two  approximations  is  equivalent  to 
assuming  that  the  unperturbed  Internal  wave  field  is  essentially 
constant  over  the  volume  of  sea  water  displaced  by  the  buoy. 
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For  spatial  wavelengths  of  principal  Interest  herein  (1,000 
to  100  meters),  this  Is  certainly  a  reasonable  assumption. 

The  technique  of  solution  involves  a  perturbation  expan¬ 
sion  In  ascending  powers  of  the  sphere  radius.  When  the  sphere 
radius  Is  small  in  relation  to  the  spatial  period  of  the  hydro- 
dynamic  flow  field,  only  the  first  term  of  the  perturbation 
expansion  needs  to  be  retained.  For  notatlonal  convenience. 

In  the  analytical  development  only  the  radius  of  the  sphere  a, 
rather  than  the  dimensionless  ratio  a/X  is  referred  to  explic¬ 
itly  as  the  small  expansion  in  parameter.  (See  discussion  on 
p.  following  Eq.  (122)). 

In  Section  II  the  general  mathematical  framework  Is  set 
forth.  The  rationale  for  the  perturbation  technique  is  devel¬ 
oped  In  Section  III  and  Appendix  A.  In  Section  IV  an  expres¬ 
sion  for  the  magnetic  field  Is  obtained  which  is  valid  in  the 
zeroth  order  approximation.  The  final  result  Is  embodied  In 
Eq.  (72),  which  appears  to  bear  a  strong  resemblance  to  the 
well-known  "cavity”  definition  of  electromagnetic  field  quanti¬ 
ties  In  an  extended  medium.  A  formula  for  the  "correction” 
term  Involving  the  Internal  wave  parameters  explicitly  is  given 
by  Eq.  (83),  while  the  corresponding  unperturbed  field  compo¬ 
nents  (l.e..  In  the  absence  of  the  buoy)  are  given  by  Eq.  (78). 
These  results  show  that  the  perturbation  of  the  field  due  to 
the  buoy  boundary  is  of  the  same  order  of  magnitude  as  the 
original  field  in  the  ocean. 

In  the  zeroth  order  approximation  the  field  is  constant 
within  the  buoy.  Thus,  In  order  to  estimate  the  spatial  gra¬ 
dients,  the  next  higher-order  term  must  be  retained  In  the  per¬ 
turbation  expansion.  The  analysis  Is  carried  out  in  Section  V. 
The  final  result  Is  given  by  Eq.  (122).  In  Section  VI  this 
formula  Is  employed  In  the  computation  of  magnetic  field  gradi¬ 
ents,  the  general  formula  for  which  is  Eq.  (124).  Despite  Its 
complexity,  the  result  admits  of  a  simple  physical  Interpreta- 
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tion:  the  effect  of  the  enclosure  is  to  symmetrize  the  unper¬ 
turbed  gradient  and  to  add  a  rotational  shear  term.  The  symme- 
trizatlon  is  in  the  form  1/2[G^^^  +  where  is  the 

gradient  of  the  q-th  component  of  the  magnetic  field  in  the 
direction  p,  in  the  absence  of  the  buoy.  As  noted  earlier, 

gg^  water,  whereas  in  the  buoy  the  gradients  must 

PQ  QP 

be  symmetric  in  the  indices  by  virtue  of  the  vanishing  of  the 
curl  of  the  magnetic  field  in  air.  The  net  effect  of  the  shear 
term  is  harder  to  Interpret.  In  some  numerical  calculations 
performed  thus  far  it  appears  to  be  small. 

Equation  (135)  translates  (12^)  into  a  formula  in  which 
the  dependence  on  the  internal  wave  parameters  is  made  explicit. 
Unlike  in  the  case  of  the  magnetic  field,  Eqs.  (78)  and  (83), 
the  formula  for  the  gradients  cannot  be  Interpreted  as  a  sum 
of  an  unperturbed  quantity  plus  a  perturbation  term.  Equation 
(135)  is  used  in  Section  VII  to  compute  the  correlation  func¬ 
tion  functions  and  spectra  for  ocean  Internal  waves.  Both 
moored  and  towed  situations  are  considered.  In  the  latter  case 
the  fast  tow  approximation  is  employed  in  conjunction  with  the 
hypothesis  of  Milder  as  given  in  [1].  The  final  formulas,  al¬ 
though  quite  cumbersome  in  appearance,  reduce  the  computation 
to  a  series  of  quadratures  involving  the  Vaisala  frequency. 

The  detailed  sequence  of  steps  required  to  Implement  the  spec¬ 
tral  calculations  on  a  computer  are  presented  in  Appendix  C. 

Only  a  few  numerical  calculations  using  the  theory  here 
developed  have  thus  far  been  carried  out .  Figure  1  shows  some 
preliminary  results  for  a  unidirectional  single  frequency  inter¬ 
nal  wave.  The  Vaisala  frequency  profile  is  assumed  exponentially 
decreasing  with  the  maximum  Vaisala  frequency  of  .833  mHz  at  the 
ocean  surface  and  a  decay  constant  of  1300  m.  A  first  mode 
Internal  wave  of  1  meter  maximum  amplitude  and  a  wavelength  of 
200  meters  (frequency  .716  mHz)  is  assumed.  The  amplitude  of 
the  Internal  wave  displacement  as  a  function  of  depth  is  indi¬ 
cated  by  the  broken  curve  (peak  at  about  100  m  depth). 
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Internal  Wave  Has  1.0-Meter  Displacement  at  Mode 


The  spectral  peak  of  the  gradient  is  plotted  for  an  integration 
time  of  1000  sec  (Vaisala  period  «  1200  sec).  The  geomagnetic 
field  is  assumed  purely  horizontal  (equatorial  zone)  and  along 
the  X  direction  (Pig.  2).  The  direction  of  travel  of  the  inter¬ 
nal  wave  relative  to  the  geomagnetic  field  is  denoted  by  w  in 
the  figure.  (w  *  0,  wave  direction  along  the  geomagnetic  field, 
w  *  Tr/2,  wave  direction  normal  to  the  geomagnetic  field).  The 
gradients  are  defined  as  follows:  component  xy  is  the  derivative 
of  the  y  (vertical)  component  of  the  field  along  the  x  direc¬ 
tion.  (The  coordinate  system  is  shown  in  Pig.  2.) 
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II.  FORMULATION  OF  THE  PROBLEM 


The  sea  surface  is  assumed  to  be  plane  and  sea  water 
occupies  the  region  y  <  D,  shown  in  Pig.  2.  At  a  depth  D  we 
assume  a  spherical  enclosure  of  radius  a,  which,  to  the  first 
approximation,  constitutes  our  model  for  a  buoy  containing  the 
magnetic  sensing  instrument  (e.g.,  a  superconducting  gradiom- 
eter).  For  the  purpose  of  the  following  analysis  the  interior 
of  the  buoy  is  filled  with  air.  Our  objective  is  to  obtain,  in 
the  interior  of  the  spherical  enclosure,  an  approximation  to 
the  magnetic  field  and  its  spatial  gradients  that  are  induced 
by  the  interaction  of  ocean  Internal  waves  with  the  geomagnetic 
field.  Analytical  results  for  computing  the  Induced  magnetic 
field  in  the  ocean  in  the  absence  of  an  enclosure  are  presented 
in  [1].  The  analysis  presented  herein  relies  heavily  on  the 
results  and  notation  employed  in  [1].  As  in  [1],  we  shall  be 
Interested  only  in  Internal  waves  with  periods  much  shorter 
than  the  inertial  period.  Consequently,  the  vortlclty  vector 
w(r,t)  can  be  assumed  parallel  to  the  ocean  surface.  We  denote 
by  V(r,t)  the  fluid  velocity  so  that 

<D(r,t)  «  V  X  V(r,t)  .  (1) 

We  make  the  usual  incompressibility  assumption  and  express  the 
fluid  velocity  in  terms  of  the  vector  stream  function  i{)(r,t), 
viz. , 


V(r,t)  -  7  X  i|)(r,t).  (2) 

•V  ^  tM 

Unless  germaln  to  the  discussion,  we  shall  henceforth  omit  the 
explicit  dependence  on  time  in  the  arguments  of  the  velocity, 
the  vortlclty  and  the  stream  function. 


n 


z 


It-t-Tt-K 

FIGURE  2,  Coordinate  System. 

In  the  course  of  the  analysis  we  shall  use  Green’s  functions 
of  the  Laplace  equation  pertaining  to  the  regions  within  and 
outside  the  sphere.  We  adopt  the  following  notation.  The 
symbols  S  and  P  will  be  employed  to  denote  the  surface  of  the 
sphere  and  ocean,  respectively.  The  region  Inside  the  sphere 
Is  designated  by  Vg.  The  water  region  between  the  surfaces  S 
and  P  Is  designated  by  Vpg  and  the  entire  undersea  region  below 
the  ocean  surface  by  Vp.  Volume  Integrals  are  Indicated  by  the 
Integration  region  Vg,  Vpg  or  Vp  and  by  the  volume  element  by 
dv.  Surface  Integrals  are  indicated,  similarly,  by  the  region 
of  Integration  S  or  P  and  by  the  surface  element  ds.  The  sur¬ 
face  of  the  unit  sphere  Is  labeled  (2  and  the  corresponding  sur¬ 
face  element  dG. 
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The  free  space  Green's  function  Is  designated  by  G^;  l.e.. 


(3) 


The  Dlrlchlet  Green's  function  that  vanishes  on  both  the  sur¬ 
faces  P  and  S  is  designated  by  G^.  The  Dlrlchlet  Green's 
functions  that  vanish  on  P  or  on  S  alone  are  designated  by  G^p 
or  Gfpg.  Similar  designations,  with  N  replacing  D,  are  used  for 
the  Neumann  Green's  function  whose  normal  derivatives  vanish  on 
the  specified  surface  or  surfaces. 


The  vertical  displacement  of  the  ocean  surface  plays  a 
very  minor  role  In  the  Internal  wave  motion  In  the  ocean.  Con¬ 
sequently,  the  boundary  condition  at  the  ocean  surface  will  be 
taken  as  V  =  0.  Also,  the  buoy  will  be  assumed  rigid  so  that 

O’ 

the  normal  component  of  fluid  velocity  on  the  surface  of  the 
sphere  must  be  zero. 


We  first  review  briefly  the  formulation  for  the  induced 
magnetic  field  In  the  ocean  in  the  absence  of  a  buoy.  The 
geomagnetic  field  is  assumed  constant  and  Is  denoted  by  As 

shown  in  [1],  the  electric  current  J,  Induced  by  the  motion  of 
sea  water,  is  given  by 


Jp  -  a  C-V<J)'  +  (B^  •  V)j|;p],  (l|) 

and  the  electrostatic  potential  in  the  water  is  expressed  by 

♦  “  -  Bq  •  ijjp  .  (5) 

The  scalar  potential  function  (^'  satisfies 

»  0.  (6) 

The  stream  function  ipp  is  obtained  from  the  vortlclty  u,  which 
we  regard  as  a  prescribed  source  function.  Since 

u  ■  7x7x^p  ■  VV»Jp  -  V^jjip,  (7) 
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we  may  use  the  gauge  condition 


7  • 


(8) 


to  obtain 


=  -  w  .  (9) 

We  solve  (9)  subject  to  the  boundary  condition  on  \|)p  that  ensures 
Vy  =  0  at  the  ocean  surface  and  is,  at  the  same  time,  compatible 
with  the  gauge  condition  expressed  by  (8).  The  boundary  con¬ 
dition  on  the  scalar  potential  <|>'  in  (6)  is  obtained  from  the 
requirement  that  the  normal  component  of  electric  current  (4) 
vanishes  at  the  ocean-air  interface.  When  the  vortlclty  u  has 
no  vertical  component,  i|)p  is  also  purely  horizontal  and  one 
finds  that 


•||-  *  0  on  P  .  (10) 

One  then  concludes  that  7^’  s  0  everywhere  and  the  electric 
current  becomes 

Jp  »  o(Bg  •  7)J,p  (11) 


In  this  special  case  the  boundary  condition  jjip  *  0  at  P  ensures 
that  V  »  0  at  P.  The  solution  of  (9)  is  then  given  by 

V 


u>  dv  . 


(12) 


One  can  show  directly  from  (12)  that  the  gauge  condition  (8)  is 
satisfied.  This  ensures  that  the  vortlclty  u  *  7x7xt|/p  with 
t|>p  as  computed  from  (12)  is  Identical  with  that  prescribed  in 
the  Integrand.  The  Induced  magnetic  field  Bp  is  determined 
with  the  aid  of  the  vector  potential  Ap,  which  satisfies 


-“o  £p> 


(13) 
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subject  to  the  gauge  condition 


V  •  Ap  »  0. 

The  solution  of  (13)  is 


(14) 


-P 


Vp*'' 
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and  the  magnetic  field  is  found  from  the  formula 

Bp  *  7  X  Ap  . 


(15) 


(16) 


We  now  consider  the  problem  of  finding  the  magnetic  field 
in  the  spherical  enclosure.  We  have 


with 


B  =  V  X  A 


A 


J  dv, 
o  -  * 


(17) 


(18) 


where  the  current  J  is  to  be  determined  from 


J  .  a[-7(>*  +  (B^  •  V)i()]. 


(19) 


The  scalar  potential  is  no  longer  zero,  since  in  addition  to 
the  boundary  condition  (10),  one  must  have 


■  ir  ’  ^?o  * 

on  the  surface  of  the  sphere.  (Here  1^,  is  the  radial  unit 
vector).  We  again  consider  only  a  purely  horizontal  vortlcity 
function.  Consequently,  at  the  ocean  surface 

i|>  ■  0  at  P. 

On  the  spherical  surface  we  require  1^  • 

i.  •  7  X  \|»  ■  0  on  S. 
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(21) 

V  ■  0,  or,  equivalently, 

(22) 


We  prescribe  the  vorticity  function  w  in  Vpg  and  determine 
from 

7  X  V  X  =  0),  (23) 

subject  to  boundary  conditions  (21)  and  (22).  An  exact  solution 
of  this  problem  is  rather  difficult.  We  shall  be  interested 
only  in  a  small  spherical  enclosure  in  which  case  one  may  employ 
the  following  approximate  solution: 


(2i*) 


Equation  (24)  gives  a  stream  function  which  vanishes  not 
only  at  the  ocean  surface  but  also  on  the  surface  of  the  sphere. 
Although  the  latter  condition  is  compatible  with  (22),  the 
divergence  of  ^  turns  out  not  to  be  equal  to  zero.  Consequently, 
the  vorticity  as  computed  with  the  aid  of  (24)  does  not  agree 
with  the  vorticity  originally  prescribed  in  the  problem.  Instead, 
one  obtains 


As  the  radius  of  the  sphere  tends  to  zero,  G^p,  so  that 

the  last  number  in  (25)  must  be  0(a).  Consequently,  for  a 
small*  sphere,  we  have 

Sactual  ”  "  ■  '26) 


I.e.,  small  in  the  sense  that  the  velocity  field  (e.g,,  the 
stream  function)  Is  constant  over  the  sphere  volume. 


We  shall  be  interested  in  computing  the  magnetic  field  and 
its  spatial  gradients  inside  the  spherical  buoy  when  the  radius 
a  is  small.  In  this  case  it  is  more  convenient  to  deal  with 
the  perturbed  quantities  B  -  Bp,  A  -  Ap,  -  \pp.  The  first  two 
are  related  by 

B  -  Bp  »  7  X  (A  -  Ap).  (27) 

Because  of  (15)  and  (l8) 


(J-Jp)G^dv  - 

~  ..r  O 


JpGodv, 


(28) 


where  the  observation  point  r  of  the  Green's  function  is  inside 
the  sphere,  i.e.,  |r|  »  r  <a. 
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III.  PERTURBATION  SOLUTION  UNDER  THE 
ASSUMPTION  OF  A  SMALL  SPHERE  RADIUS 


The  free  space  Green's  function  given  by  (3)  has  the 


power  series  expansion,  when  r  <  r', 

n«0 

where  is  the  Legendre  polynomial  of  order  n  and  y  is  the 
angle  between  the  vectors  r  and  r’.  It  follows  that  for  the 
first  term  on  the  right  in  (28) 


(29) 


'‘c  /  <£- 


Where 


'PS 


n*0 


(30) 


„  f  J-Jp 

-n  *  w  y  fnrfT  • 


(31) 


'PS 


In  (30)  and  (31)  dv'  has  been  written  in  place  of  dv  to  empha¬ 
size  that  integration  is  with  respect  to  primed  variables  asso¬ 
ciated  with  G^  and  that  the  result  of  the  Integration  is  still 
a  function  of  the  unprimed  variables.  Since  P^  ®  1  it  should 

be  noted  that  a  =  constant.  It  follows  that 
~o 

(J-Jp)  Gq  (Vx  (32; 

n*l 

In  accounting  for  the  second  term  on  the  right  of  (28)  the 
identity 


''ps 


(33) 


15 


implies  that 


/jp  G  dV  =/ 7'xJpO  dv'  -  ft 
Vo  °  Vo  °  s 


X  JpGods'. 


Then,  using  the  power  series  expansion  (29)  for  when  r  <  r' 
and  the  expansion  obtained  by  interchanging  r  and  r'  when  r  >  r', 
it  follows  that 

^^/*£p  ^  Y)dr*dJ2 

V  Q  r 


dr 'dS2 


Q  0 


00  ^ 

-  (l)  f  <lr  ='  ip>lr-a 

n=0  n 

“  m  f  Ij.,  X  Jp(0)P^(cos  Y)d:i+0(a^),(25) 

n=0  '  '  fi 

2 

where^O(a  )  has  been  written  in  place  of  functions  of  the  form 

a^  X  >  C  (^)^.  In  deriving  (3^)  use  has  been  made  of  the  con- 

tinulty  of  Jp  to  relate  Jp  on  S  to  its  value  at  the  center  of 
the  sphere,  l.e.. 


JpL.^a  =  £p^°^ 


It  follows  from  (27),  (28),  (32),  and  (35)  that 


H  u.- 


*  (B-B  )  »  Z-f  (Ip  •  Vxa  )r' 
~r  ~  ~r  n 


+  Uo  W  ^  (|)V  ir  •  Ir-  *  Jp(0)P„(dosY)dn  +  OCa^) 


00  OO 

.^(1^  •  ir’'ir'''n<=°=  I-)" 

n»l  n=0  Q 

+  0  (a^) 

00 

=  (1  •  7xa^)r^  +  0  (a^).  (36) 

-r  ~n 

n=l 

The  last  relation  follows  from  the  fact  that  the  integral  in 
the  penultimate  relation  satisfies 

ir  X  ij,,  P^(cos  Y)dfi  =  0. 

Q 

This  can  be  seen  from  the  fact  that,  if  i^,  is  taken  to 
polar  axis  relative  to  the  primed  coordinates  in  (37), 
rotating  coordinates  in  the  integrand,  then 

1  X  1  ’  =  i . ,  is  the  azimuthal  unit  vector 
^r  ~r 


(37) 

be  the 
l.e.,  by 


and 


P^(cos  y)  =  P^(cos  0).  Since  dfi  =  sined0d(t>  one 
readily  verifies  that  integration  over  a  complete  sphere  yields 
zero. 


According  to  (36),  on  the  spherical  surface. 


1  •  B 

ir  ~'r»a 


(ir  •  ^  ^  “n^r-a^""  ^  ^ 

n*l 


*  b^(0,(Ji)  +  0(a  ). 

Since  inside  the  sphere 

V  X  B  -  0, 

there  exists  a  scalar  function  ♦  in  terms  of  which 


(38) 


B  •  V  9 


(39) 
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inside  the  sphere.  Also 

V  •  B  *  7^  $  =  o  (40) 

inside  the  sphere,  and  (38)  and  (39)  imply  that  on  the  surface 
of  the  sphere 


To  find  the  magnetic  field  within  the  spherical  enclosure, 
it  is  only  necessary  to  solve  Laplace’s  equation  for  $  in 
subject  to  the  boundary  condition  (4l)  on  S.  The  solution  for 
small  a  does  not  vanish  in  the  limit  as  a  approaches  zero,  but 
satisfies  the  Laplace  equation  boundary  value  problem  with  the 
radial  derivative  of  $  equal  to  b^(6,(|)),  a  quantity  given  by 


b  (0,(|))  =  lim 
a-*-0 


ir  *  ?P'r 


=a  *  E  <ir  • 


n=l 


(42) 


In  calculating  the  magnetic  field  inside  the  sphere  for 
small  a,  it  is  appropriate  to  use  a  power  series  in  r,  or  rather. 
Accordingly,  the  potential  function  i  is 

cL 

*  ■  E  HI) 

n=l 


where  the  Z^(<|),0)  are  the  terms  in  the  spherical  harmonic  ex 
pension*  of  the  boundary  value: 


The  Z  (<>,0)  denotes  the  sum  d^P^(cos0)e^'"‘**,  where  d  are 

n  m  n  *  m 

m“— n 

suitable  constants  and  P""  (cos0)  are  associated  Legendre  poly¬ 
nomials  .  “ 
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94> 

3r 


n=l 


b^(^,e) 


(4A) 


The  fact  that  the  series  in  (43)  and  (44)  are  both  missing 
the  term  corresponding  to  n=0  Is  a  consequence  of  Gauss'  theorem 
and  the  fact  that  the  divergence  of  any  magnetic  field  is  zero. 
That  iSj  the  Integral  over  S  of  the  magnetic  field's  normal  com¬ 
ponent  vanishes,  but  the  same  integral  provides  the  coefficient 
of  the  zeroth  order  term  in  the  normal  component's  spherical 
harmonic  expansion. 

According  to  (43), 

B  -  7  4 . 2  i  7jj  zj  (l)"*''  .  iH5) 

n=l 


In  (45)  the  operator  is  the  angular  part  of  the  gradient 
and  is  defined  by 

7  si  ^  i _  i 

Q  ie  ae  sine  i(})  3(fr 
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IV.  EXPLICIT  EVALUATION  OF  THE  MAGNETIC  FIELD  TO  THE  ZEROTH 
ORDER  IN  a 

The  coefficients  In  the  power  series  (45)  for  the  magnetic 
field  are  obtained  from  terms  In  the  spherical  harmonic  expan¬ 
sion  (44)  of  the  function  b„(0,4)),  given  by 

3L 


b^C^.e)  -  1^  •  (1^  • 


('16) 


n=l 


The  first  term  on  the  rleht  of  (46)  Is  oresumably  known,  and 
Is,  In  fact,  to  the  lowest  order  In  a,  the  radial  component  of 
the  unperturbed  magnetic  field  at  the  center  of  the  sphere. 

The  remaining  terms  must  be  obtained  from  the  vectors  defined 
by  (31). 

According  to  (19)  and  (11),  the  current  appearing  In  (31) 

Is  given  by 


(47) 


where 


a  (B^  •  V)  (Vf-jt'p) 


and 


Jg  «  -  a  V  . 

It  will  now  be  shown  that  the  term  contributes  a  term  of 
order  a  to  the  magnetic  field  Inside  the  spherical  enclosure. 
It  will  be  found  later  that  the  contribution  of  Is  of  order 
one  In  a. 
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Because  of  (27),  (28)  and  (29)  the  contribution  of  Jg  to 
the  magnetic  field  is 


?E 


-yoC  7 


1 


V(|)’ 


G  dv' 
o 


VG  xV(t)'dv* 
o 


=  f  7*G^xV(J)’dv’  =  y^o  f  V' x(G^V(|)' )dv’ 

*'V  •At 

'^ps  ''ps 

=  f  n*xV(J>*G^ds*  . 

•'p+s 


(^♦8) 


The  function  <|)'  in  (48)  is  the  solution  of  the  boundary 
value  problem  given  by  (6),  (10),  and  (20)  in  which  the  normal 
derivative  of  <|»*  is  given  on  S  and  P.  It  is  shown  in  Appendix 
A  that  the  solution  of  a  similar  problem,  in  which  the  function 
rather  than  its  normal  derivative  is  prescribed  on  S  and  P  can 
be  obtained  approximately  by  ignoring  the  boundary  condition 
on  P.  The  approximate  solution  is  correct  to  the  lowest  order 
in  the  ratio  A  similar  argument  can  be  given  for  the  case 
in  which  the  normal  derivative  is  prescribed  rather  than  the 
function.  Accordingly,  it  can  be  asserted  that  the  boundary 
value  problem  determined  by  (6)  and  (20)  gives  ({>'  correctly 
except  for  terms  of  order 

Thus,  ())'  can  be  written  in  the  form 


4,* 


.’f' 


Zn(*,e)+0(g), 


W) 
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where  the  Z^((|),9)  are  the  terms  in  the  spherical  harmonic  ex¬ 
pansion  of  the  boundary  value  on  S: 


li’  I 
3r  'r=a 


i 


~r 


i 'r»a 


Then  V())*  has  the  form 

»  ^  /_»n+2 

~2-(f)  in'*-®'*  “  §). 

n»0 

where 


(50) 


5„  (♦.6)  -  Z„(«.e)  ir  -  ^  Vn'*-«>- 


Neglecting  terms  of  order  ^  in  (50)  is  equivalent  to  drop¬ 
ping  the  Integral  over  P  in  (48).  Then 


?E 


X  V(|)’G^ds». 


(51) 


Prom  (50),  (51),  and  the  power  series  expansion  (29)  of  G 
for  an  observation  point  Inside  the  sphere  it  follows  that 


(cos  y)  1^,x  Cj^((J)' ,6' )dn. 


(52) 


Since  in  (52)  r  <  a,  it  follows  that  Bg  is  at  least  0(a),  as 
was  to  be  demonstrated.  Consequently,  except  for  quantities 
that  contribute  terms  of  order  a  to  the  magnetic  field 


J  -  Jp  ~  a  (B^  •  (53) 

We  must  now  examine  the  representation  of  the  stream  func¬ 
tions  J  and  )|»p.  Since  G^  =  Q^p  +  Gp  -  G^p,  Eq.  (24)  may  be 
written  as  follows: 

■  /  Gpp  u  dv  ^  AQudv,  (54) 

Vps  ' 
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where  AG  is  the  difference  between  the  Dlrlchlet  Green's  func¬ 
tion  for  the  region  with  boundaries  P  and  S  and  the  Dlrlchlet 
Green's  function  Gj^p  for  the  region  with  the  boundary  P  alone. 
Similarly,  we  write  the  integral  (12)  as  a  sum  of  two  parts: 


/  OdpS’*''  -  f 


(55) 


'p  ’PS 

Subtracting  this  from  (5^)  we  have 


\|»  -  /  A  G  03  dv  -y 

~  “  Vps  ~  *^S 


Gjjpudv . 


(56) 


It  is  convenient  to  express  the  various  Green's  functions 
in  terms  of  a  radius  vector 

r(x,y,z)  =  xi^  +  yiy+^iz 
and  Its  Image  r  In  P, 

r  =  r  (x,2D-y,z). 

Clearly,  on  P 


r|p  =  r(x,D,z)  *  r(x,D,z)  =  rip  (57) 

Thus, 

°DP  ■  0o'r>r'>-0o<r>j:’>>  '^e) 

since  (57)  Implies  that  G^p  defined  by  (58)  vanishes  on  P,  both 
terms  on  the  right  of  (58)  satisfy  Laplace's  equation,  and  G^p 
has  the  same  singularity  as  G^  when  r  approaches  r'  in  Vp. 

Making  use  of  the  expansion  (29),  one  can  obtain  an  esti¬ 
mate  of  the  second  term  on  the  right  side  of  ( 56)  by  observing 
that 


/*Gppa)dv’»  Y)dV-  ur'"p^(cosY)dv'l 

Vo  n»0  ^  Vo  ^  ir  -* 


=  t?  £  H^(f)  / Y)dn 


-(f  Y. .. 


(cos  $)dG  +  O(a^) 


+  0  (a^). 


where  w  Is  the  value  of  u  at  the  center  of  the  sphere  and  where 
““0  ~  ^ 

Y  is  the  angle  between  r  and  r'  and  ^  is  the  angle  between  r  and 

r’.  Because  of  (59)j  the  second  term  on  the  right  side  of  (56) 

"  2 

will  be  carried,  simply,  as  a  term  of  order  a  in  comparison 

with  the  first  term  whose  order  will  also  be  obtained. 

It  is  shown  in  Appendix  A  that,  to  the  lowest  order  in  sr. 


^a  -2a  rr'cosY-tr  r’  ^a^-2a^rr’cos9+r  r’ 


F==l- 

rr’cos^+r^r’ 


(60) 


Since 


and  in  V, 


r  >  a 


AO  can  be  expanded  in  the  power  series 
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Hence, 


At 

P^(cos  Y)dv 


(cos  y)dv 


dv*  +  O(a^) 


Hence,  to  the  lowest  order  In  a, 

^  -  ij/p  ~  a  (^  -  p)iM, 

where  y  is  the  constant  vector  given  by 


According  to  (53)  and  (62) 


(62) 


(63) 


?ojif\y,  (6i|) 

where,  as  usual,  the  vectors  i^  and  are  unit  vectors  in  the 
directions  indicated  by  their  subscripts.  Substituting  (6i|)  in 
(31)  yields 


go  •  V 

r’2  r.n+iy  n 


(cos  Y)dv*  y. 


C65) 


Since  terms  of  order  ^  are  to  be  neglected,  the  integration 
region  V„<j  can  be  replaced  by  the  region  between  S  and  a  sphere 
With  the  same  center  and  radius  D.  Then  the  term  involving  r 
can  be  neglected  as  well.  The  result  will  be 


-  ^7/ 

n  a 


!o 


:;n+r 


P^(cos 


Y)dr 'dfly. 
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(66) 


which,  after  neglecting  terms  of  order  ^  becomes  for*  n  /  0 
a  ~  _ _  f  H- 

The  integral  on  the  right  side  of  (66)  can  be  simplified 
by  rotating  coordinates,  if  necessary,  so  that  the  polar  axis 
of  the  coordinate  system  corresponding  to  the  integration  vari¬ 
ables  has  the  direction  1^.  Then  the  angle  y  becomes  0*.  If 
we  express  the  geomagnetic  field  in  the  form  =  1^  + 

^B^  -  B^  [cos  6  1^  -  sin  6  1a1»  the  quantity  B^  •  i^t  may  be 
o  or  ~o  ~<p  '  ^  ~r 

written  as  follows: 


?o  •  i  V^o^-®or  (<!>’+ 6 )+B^r°°s  0’ 


Thus,  (66)  becomes 


a  — 
~n 


<Tu  f2ir  /‘n/r 

{(^ 


Bo^-Bor^  sin  0’  cos(<|)’+B) 


ov 

n2a 


cos  0'1p  (cos  0’)sln  0*d0*d^’ 
or  in 


(cos  0')  sin  0’  cos  0'  d0’ 


ou 

n2a 


-1 


)dx 


1.  e. , 


^^o^or  u;  n  »  1> 


3 

0 


J  n  +  1, 


?1 - 3"  ^ir  *  ?o^y*  ^^7) 

5n  ~  0,  n  +  1. 

•a.  does  not  enter  into  the  calculation  of  the  magnetic  field. 
See  (32) 
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Then,  according  to  (67) » 


7  X  ~  V(i^  • 

,n  +  l. 


But 

’'So'lr'  ■  r  ’'So  •  ^  '?o  ■ 

■  ?  <?o-®orir)- 

Prom  (68)  and  (69)  it  follows  that 

ay 

ir  *  ^  3l~  IF  ir  ""  ®o  *  y  » 


(69) 


(70) 


Ir  .  7  X  On  ~  0,  n  +  1- 


Thus,  for  the  boundary  condition  determined  by  (38),  it 
follows  from  (70)  and  (4l)  that,  to  the  lowest  order  in  a. 


3«i 


ba(4',0)~^lr.a  ~  ir  *  ?P'r-a  *  ~T  -r 


1_  X  •  V' 


av 

~  1  •  B„(0)  +  1_  X  B. 


tr  ~P’ 


“T  ir  ^  eo 


-  1^  •CBp(O)  +  -3^  X  y],  (71) 

where  Bp(0)  Is  the  unperturbed  magnetic  field  at  the  center  of 
the  sphere.  According  to  (71)  the  normal  component  of  B  on  S, 
to  the  lowest  order  In  a.  Is  equal  to  the  normal  component 
•  P  of  a  constant  vector  S  defined  by 

X*  ^ 

ay 

6  -  Bp(0)  +  -5^  X  V  . 


) 
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The  vector  B  whose  divergence  and  curl  vanishes  In  S  and 
whose  normal  component  on  S  is  1^  •  B  is  the  constant  B,  Itself. 
This  can  be  verified  algebraically  by  expanding  bg^((^,0),  given 
by  (71),  in  a  spherical  harmonic  series  and  substituting  into 
(45).  In  this  case  the  series  consists  of  a  single  term 
Zj  (<(>,6),  the  result  of  calculating  the  radial  component  of  a 
constant  vector.  Thus,  to  the  lowest  order  In  a,  we  have  the 
final  result 

B  -  Bp(0)  +  —  B^  X  u,  (72) 

where  ' 


Equation  (72)  states  that  the  magnetic  field  Inside  a 
spherical  buoy  of  small  radius  Is  spatially  Invariant  and  is 
given  by  the  sum  of  the  unperturbed  field  (as  computed  at  the 
origin  of  the  coordinate  system  of  the  sphere)  and  a  correction 
term  which  Is  proportional  to  a  volume  Integral  comprising  the 
vorticlty  function.  The  fact  that  the  field  Inside  the  spherical 
cavity  Is  spatially  Invariant  Is  analogous  to  the  situation  that 
arises  when  a  dielectric  sphere  Is  placed  In  an  electrostatic 
field:  there  the  field  Inside  the  sphere  is  also  a  constant. 
Indeed,  the  presence  of  the  factor  1/3  In  the  correction  term 
Is  reminiscent  of  the  depolarization  factor  for  the  electro¬ 
static  field  measured  In  a  hollow  spherical  cavity  in  an  ex¬ 
tended  dielectric  medium.  Thus,  when  n  Is  the  (relative) 
medium  dielectric  constant,  the  electrostatic  field  within  a 
spherical  cavity  Is 

where  Is  the  Incident  (unperturbed)  electrostatic  field. 
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We  now  relate  Equation  (72)  to  the  parameters  of  the  in¬ 
ternal  wave  field.  We  denote  the  cartesian  coordinates  with 
origin  at  the  ocean  surface  by  where  n  ^  0  Is  occupied 

by  sea  water.  These  coordinate  axes  are  chosen  parallel  to  the 
x,y,z  coordinate  employed  In  the  preceding  discussion  (see  Fig. 
2).  We  denote  the  Instantaneous  vertical  displacement  due  to 
internal  wave  motion  by  q(C,n,?,t).  The  horizontal  coordinate 
pair  (C>?)  will  also  be  denoted  by  the  vector  ic, 

<  «  C  +  C  • 

We  now  Introduce  the  two-dimensional  Fourier  transform  represen¬ 
tation  of  q: 


00 

q(!c,n,t)  =  y*  J  q(K,n,t)d^K,  (73) 


where 


^  nC,(K)  a;;(k)  -in„tl 

q(K,n.t)  e 

L  n  In  J 

The  4>y^(h)  are  internal  wave  mode  eigenfunctions  satisfying 

^  <fr„(Ti)  +  K^— -  lVr,(b)  -  0, 


(7^) 


(75) 


dn 


n 


with  the  normalization 


y  <|>j^(Ti)N^(Ti)dn  •  .  (76) 

A  typical  number  of  the  sum  In  (7^)  can  be  Interpreted  as  a 
simple  harmonic  Internal  wave  with  angular  frequency  travel¬ 
ing  in  direction  K/K,  as  Indicated  by  the  angle  w  in  Fig.  2. 
Each  component  Internal  wave  Induces  a  corresponding  partial 
magnetic  field  waveform  which  we  denote  by  ^(K,n,t).  The 
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actual  field  is  then  given  by  the  superposition 

B(n,n,t)  f  ^  B^"^(K,n,t)d^K  .  (77) 

In  the  following,  we  shall  be  dealing  explicitly  only  with  a 
typical  member  and  use  the  notation  Bp”^  for  the  unper¬ 

turbed  field  in  the  ocean.  Prom  [1],  p.  5^1,  one  has 


n(n)  _  p(n)  .  g(n) 
?P  ■  ?PH  ^  U  ®Pn 


?PH  \ K 


?0  0„(K.t). 


Where 


'l: 

1 


^  \  ^  y  f  mm 

kI 

ry^n  *  "vf/tj-  ?0 


.  ifl^t 

U  (K,t)  -  (K)  e  ^  +  a;  (K)  e  ^  , 

n  n  -•  n  ^ 


"  •/n 


^n  “  ® 


~  ♦nCnMdn’. 


(78a) 

(78b) 

(78c) 

(78d) 

(78e) 

(78f) 


The  subscript  H  in  (78a)  and  (78b)  is  used  to  denote  the  fact 
that  this  vector  is  purely  horizontal.  These  expressions  give 
the  unperturbed  magnetic  field  below  the  ocean  surface.  We 
now  obtain  the  perturbation  due  to  the  spherical  enclosure,  as 
given  by  the  second  term  on  the  right  of  (72).  We  first  com¬ 
pute  y.  When  the  location  of  the  center  of  the  sphere  is  ex- 
pressed  in  terms  of  5,ti,C  (n  ■  -D  <  0)  the  Integral  may  be 
written  as  follows: 
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'  ■  ^  r 


fe)(x+g.  y-Hi.  Z+C) 

dxd z 

x'^  +  y'^  +  z 


^  r  ^^'Cr 

y-oo  J-^J.co  y(n'-n)  +(C’- 


O^+Cc'-c)^ 


Using  the  identity 
1 


/  ='■■" V  ■  - -  _J^  /‘“/’“d^K  _^-K|n-n' 

Tty  (n-n')^  +  U-Vr  +  (2it)  J_„J.o.  2k 


-iK 


we  obtain 


if  =y*y*  d^K  ^  ~y*°  dn’^Cn’.K) 


where 


u(n*>K)  =  — ^-g-  oaCn’ ,<»  )d^K’ 

(2Tr)‘^  y.cojLco 


Prom  [1]  one  finds 


_  2 

u(ti’,K)  =  i(l^  xK)^^  N_(|1I  <|»^(n’ )U^(K,t) 


TT  «n 


=  2^  u^"^n’,K). 
n 


With  the  aid  of  (8l),  (79)  may  be  written  as  follows: 


!  ■££  A  L‘  !. 

n 


•  (k-k'  ) 


(79) 


(80) 


(81) 
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»».  -  ... 


f 


i 


> 


I 


? 


where 


1  X  K  /u^(K.t)\  I  . 

i„  -  1  e-''!"-''  I  N2(n')*„(nMdn’. 


Prom  the  differential  equation  (75)  it  may  be  shown  that 

'  N2(n' 


Kn  (82) 


Thus,  the  correction  term  that  must  be  added  to  (78)  to  account 
for  the  effect  of  the  enclosure  is  given  by 


?o  *  !in  ' 


lou. 


T” 


5  •  B  \  1  -  (i  •  B  )  5 

K  ‘'®)  T 


^<l>n(n)  4'^(0)e 


Kn, 


K 


2K‘ 


1.  (83) 


Comparing  this  expression  with  (78)  we  observe  that, 
generally,  the  "correction”  term  is  of  the  same  order  of  magni 
tude  as  the  unperturbed  field. 
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V.  THE  MAGNETIC  FIELD  TO  THE  FIRST  ORDER  IN  a 

Since  In  the  zeroth  order  approximation  the  magnetic  field 
in  the  spherical  buoy  Is  constant,  the  (spatial)  gradient  of 
the  magnetic  field  is  zero  in  this  approximation.  To  obtain  a 
numerical  estimate  of  the  gradient,  the  magnetic  field  must  be 
computed  to  the  first  order  in  a. 

The  magnetic  field  Is  determined  by  the  currents  and  ^ 

as  indicated  in  (^7).  It  is  shown  in  Appendix  B  that  the  com- 

blnation  of  Jg  is  0(a).  Consequently,  to  find  the  field  to 

within  0(a)  we  need  consider  only  The  magnetic  field  is 

determined  by  the  a„  defined  by  (31),  since  it  has  already  been 

~n  2 

shown  that  the  integral  over  Vg  contributes  a  term  of  order  a  . 

According  to  (53),  (56),  and  (61)  the  current  that  determines 

the  depends  upon  vectors  defined  by 


,n+l 


P^(cos  Y)dv ' . 


(84) 


The  vectors  H  can  be  written  in  the  form 
~n 


» y  p^(cos  y)  p^((^»,e’)dn. 


where 


D  sec  0 


(86) 


It  is  clear  from  (86)  that  for  n  =  0  or  1  is  0(1)  in  a,  but 

for  n  =  2  F  is  0(log  a).  For  larger  values  of  n,  i.e.,  n  >  2, 

F  is  0  (a^"*^^. 

~n 

Each  cartesian  component  of  F^^  can  be  expanded  in  a  sur¬ 
face  spherical  harmonic  series.  The  results  leads  to  a  vector 
expansion 


GQ 


where  is  a  vector  each  of  whose  cartesian  components  is  a 
spherical  harmonic  of  degree  m.  It  follows  from  (85)  and  well 
known  properties  of  surface  spherical  harmonics  (cf.  MacRobert, 
Ch.  VII)  that 

5n  ■  Mn  ^88) 

Where  has  been  written  in  place  of  Z 

-n  ^  ~nn 

From  (56)  and  (62),  neglecting  terms  of  order  a  and  of 
order  it  follows  that 


CO  n  / 


From  ( 53)  it  follows  that 


J  -  Jp  ~ 


^2n+l 


where 


)Z„<*,9)L  (91) 

is  a  vector  each  of  whose  cartesian  components  is  a  surface 
spherical  harmonic  of  degree  n+1.  This  last  statement  follows 
from  the  fact  that  each  component  of  each  term  in  the  series 
(89)  for  f ~ separately  satisfies  Laplace's  equation  and 
still  does  so  after  the  operator  (B„  •  V)  is  applied. 

•K  0 


According  to  (31),  then, 


2n 


v-^ 


Utt 


m=o 


2m+l  z'  ~m+l 


V 

^PS 


v„x-,U'.e') 


,,m+n+3  n 


P  (cos  y)dv' 


ay  -ss_  o  XI  /•  Z*^  ® 

■  ThT  2-  ^  ^7  dr'  an 

m=o  4  _,m+n+l 


a  ^,m+n+l 
r ' 


a^  ^  / 

TT  m+n  J 


Y)dn 


m=o 


n 


-  "^^0  Y  (<(.,0),  (92) 

—  ^  -n 

4n  -1 

wherein  terms  of  order  ^  have  been  dropped  and  the  same  prop¬ 
erties  of  spherical  harmonics  used  in  deriving  ( 88)  have  been 
applied.  It  should  be  noted  that  for  n  ®  1  or  2  is  0(1)  in 

a,  for  n  =  3  Y  is  0  (log  a),  and  for  n  >  3  Y  is  0(a^“^). 

-^n 

According  to  (88) 

Zl(<l>.e)  =  Hi,  (93) 

and  therefore,  according  to  (84) 


Zi((|),e) 


cos  Y  dv’ 


cos  Y  dV . 


(94) 


Thus , 


(♦jQ)  ~  y^ 
where 


cos  9  +  (jii  cos  ^  +  y^  sin  4i)  sin  0, 
-3  "  ^  cos  e*  dv’. 


Vp  r 


(95) 

(96) 
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and 


sin  6'  dv' , 


(97) 


Then 


^2 


sin  4i '  sin  0 *  dv' . 


(98) 


B  )cos  e  +  (i  •  B  )sin  sjv- 
—  o 

B^)  cos  (f>  sin  0  -  (io  •  B  )cos  (J»  cos  0 

w  -vO 


+  9  I  Ut  + 


j  +  1^2 (ij,*  B^)sin  (j)  sin  0 
-  (io*  B„)  sin  (J)  cos  0  -  (i.  *6^)008  (()  v-* 

•*0  -0  -(p  -0  J 


(99) 


According  to  (38),  if  the  boundary  term  b_((fr,0)  is  supple- 

a 

mented  by  a  quantity  A0  given  by 


A6 


V  X  012) 


r=a 


(100) 


because  of  (92)  and  the  order  in  a  of  Y  for  n  >  2,  the  terms 

<3  n 

neglected  in  b^(<f,e)  will  be  0  (a-^  log  a)  and  higher.  In  order 

cl 

to  obtain  the  complete  field  correct  to  terms  of  order  a,  how¬ 
ever,  it  is  also  necessary  to  replace  the  quantity  B_(0)  used 

••  r 

in  (71)  and  (72)  by  B„|  correct  to  the  first  order  in  a. 

JT  r*ci 

Thus,  the  new  boundary  value  b  ((J»,0)  is  given  by 

a 


3  Bp 
-r  I 


OV 


b^((>,0)~l^-Bp(O)+  3:r‘?o*  H  ^  ir*  ^Ir-O"^  Ir-a 

(101) 

+  0  (a^  log  a). 
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wherein  the  fourth  term  on  the  right  side  of  the  equation  re¬ 
sults  from  (92)  and  (100). 


The  fourth  term  on  the  right  side  of  (101)  is  actually  0(a) 

p 

rather  than  0(a)  because  of  the  curl  operation.  The  expres¬ 
sion 

r  X  V  •  [r^  y2(<))  0)]=  r^i^  •  V  x  ^102) 

is  a  homogeneous  polynomial  of  degree  two  in  the  cartesian 
variables  x,  y,  and  z  because  the  cartesian  components  of 
y2((|>,6)  are  all  surface  spherical  harmonics  of  degree  two.  It 
then  follows  (cf.  Ref.  3,  p.  1^0)  that  the  fourth  term  on  the 
right  side  of  (101)  can  be  expressed  as  the  sum  of  a  surface 
spherical  harmonic  of  degree  two  and  a  surface  spherical  har¬ 
monic  of  degree  zero,  l.e.,  a  constant.  That  constant  is  the 
Integral  over  the  unit  sphere  of  the  original  term,  for  which 
according  to  Gauss*  theorem 

.^i^-V  X  V  •  (V  X 

•'s  '  '  \ 


Hence  the  original  term  is  a  surface  spherical  harmonic  of 
degree  two. 

The  third  term  on  the  right  side  of  (101)  can  be  written 

^  ir 

\  ! 

(‘r  •  ^lr.o)(lr  ’  1.)] 

j^B2(<>,e)  +  a  (lOH) 


3Bp 

/ 

\ 

m 

3r  'r«0 

^(ir  • 

plr.O 

3Bp  \i 

f  ) 

i  ( 

1  •  i 

3x  'r®0ll 

Ur  ~xJ 

'  V~r  3y 

a 
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where  B2(<)),6)  can  be  recognized  as  a  surface  spherical  harmonic 
of  degree  two  because  of  the  Identity  (B-4).  Here,  of  course. 

Is  given  by 

3B 

B2(<|),e)  =  Ij,*  (105) 

These  considerations  make  It  possible  to  write  an  expres¬ 
sion,  correct  to  terms  of  order  a^  log  a,  for  $,  the  function 
that  satisfies  Laplace's  equation  subject  to  the  condition 
that  on  S  Its  radial  derivative  Is  equal  to  the  function  bg^((J),e) 
given  by  (101).  Since  the  magnetic  field  Inside  the  sphere 
Is  Vi,  It  Is  then  given  approximately  there  by 


B  ~  Bp(0)  + 


ay 


av^ar 

°  ?0  -  S  ^  15^ 


|(ir  *  5  ''«)  Ir  •  I 


+  r 


(ir  *  I 


r«a' 
(106) 


correct  to  terms  of  order  a. 

Because  of  (lO^l)  It  follows  that 

x(6.  •  r)  +  z(6  .  r)  +  y(B.  •  r) 

B„((|),e)  =  _ '2  -  -3  ~  ^ 


where 


3  Bp 

~1  r*0, 


^?P|  8  .  ^?P| 

3z  •r-O  ,  Z3  3y  'r»0  ’ 


(108) 


Then,  according  to  (106),  the  contribution  of  B2(4>,e)  to  B  can 
be  written 


/  1  V  /  1  \  +  z  e?  +  r 

{ir*  i  •  Kir"  J 


(109) 


The  vector  operator  In  (109)  can  be  written 


1  V  ^ 

2 


J  r  V  + 
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.  1  7  +  .  1  r  •  v)r 

7  .  i.  (l  -  I  .  .  7)]r. 


(110) 


It  is  now  useful  to  observe  that  relations  such  as 


and  thus. 


hold. 


r  (111) 


(112) 


Relations  similar  to  (111)  and  (112)  can  be  used  to  evalu¬ 
ate  (110).  Thus,  it  follows  that 


*  I  '’a)  -  l|[ixS  *  (ix-  r)  j]-?] 

*  [lx?  *  {h-  r)i]  -h 

*  [iy?  *  (iy-  r)i] 

where  I  is  the  unit  dyadic  defined  by 


I-ll  +11  +11. 

t  ix~x  -y~y  -z-z 


(113) 


(11^1) 


It  follows  from  (91)  and  (95)  that,  except  for  an  addi¬ 
tive  constant, 

Y2(*,e)  -  U^((|),e)W;L  *  U2(*,e)u2  +  U3((J),e)y3,  (115) 
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where 


U3^(<|),e)  »  I 


4  *^?o* 


,  U2(<|),e)  =  ^ 


4  ‘ 


(116) 


4  y^®o  ‘ 

Un(<|),e)  *  -F  y-  "  . 

J  ^  T 

From  (115)  it  follows  that 


(117) 


i  •VxYo*-(rxVU*y,  +rxV  U„‘  +  r  x  V  •  y,). 

~r  ~2  r  1  ~  d-d  -  5  -i 

(118) 


From  (ll6)  It  will  be  found  that 


„  „  -  n  r/So  •  r 

1  ■  5 


X  r  X 


r  X  + 

-»A 


(119) 


Similar  relations  hold,  of  course,  for  r  x  7U«  and  r  x  7U-. 

"•  t  ^  0 

It  then  follows  that 


[••  Ir  • 


»  r  1  •  V  X  Y,  . 

-r  -2 


(120) 


From  (110),  (119),  and  (120)  It  follows  that 

{ir  *  i  Vn)(ir  •  I  [I-  ^  ?  •  ’)] 
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W'  ■  -■ 


y  r  X  B 


“]■ 


(121) 


■  -(ix  >=  iz  ^  ^-2*  iy  “3)1  ?o-^  (!o- 


[(r  •  ?o  ^  Hi)lx*^(!o  *  !;i)+(!:  •  ?o^i2)iz*  ^(?o='  Ha) 

+  (r  •  !o='  H3)vK-»  -3)]} 

-  5I?{[?o(h^1x)*  ix(r*!o)*  (ix-  r)(i  >'  ?o)  *  (?o-  H)ii]-  Hi 

*  [?o  (r "  iz) "  ix(r "  ?o)  *{iz  ■  r)(i  ='  !o)"(?o  •  r)i2]  -h 

(121) 

*  [?o  (h  ==  ly)  My  (e  -  ?o)*(iy  ’  E)(i  ?o)"(!o  ’  r)i3]  -H 


where 


3 

(121a) 


;1  *  iz~y  “  iyiz*  ~2  “  iyix  "  ix~y»  ~3  *  -x~z  “  ~z~x’ 


Expressions  (113)  and  (121)  can  be  substituted  into  (106) 
to  obtain  an  explicit  expression  for  the  magnetic  field.  Then, 
with 

!o- 

the  magnetic  field  can  be  written: 

5  •  [Sb*  ^  (?o'‘  h)  *  I  |[ixE  ♦  m]  •  h  *  [izE  ♦  h}  h 
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*[iyr  *  yj]  •  Sjj  +  ^|[?o(r  *  lx)  *  lx(!:  *  ?o)*  *i*!o 


The  foregoing  analysis  shows  that  terms  neglected  In  the 
approximations  for  the  magnetic  field  and  its  gradient  are  of 
the  same  order  as  a^  multiplied  by  the  quantities  given  by 
(8^1).  That  is,  the  terms  neglected  are  of  the  order  of  a’^  times 
higher  derivatives  of  the  quantity 


G^wdv* . 
0" 


As  Indicated  In  (79),  the  convolution  theorem  implies  that  the 
horizontal  Fourier  transform  y  of  y  has  the  form 


y(K) 


K)  ^-1  K 


K  -  Kln-n* I dn 


It  follows  that  an  n^^  order  derivative  of  y  has  a  Fourier 
transform  that  Is  proportional  to  K”y(K).  Thus,  a  condition 
that  might  be  Imposed,  Justifying  the  magnetic  field  approxi¬ 
mation  used  In  this  paper,  can  be  stated  as  follows:  the 
vorticlty  spectrum  Is  negligible  except  for  the  wave  number 
region  In  which  Ka  <<  1.  Because  of  the  monotonlc  relation 
between  frequency  and  wave  number,  this  condition  Is  consistent 
with  the  static  approximation  used  In  calculating  the  magnetic 
field. 
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VI.  THE  MAGNETIC  FIELD  GRADIENT 


The  magnetic  field  gradient  relative  to  unit  vectors 
Ip  and  Iq  is  defined  by 


If  (123) 


where 


is 


(123) 


applied  to  (122)  the  result  is 

«p,  "  ?  (il  •  §1  ^  52  •  §2 

^  -75^  (El  •  ti  ^  12 


^  53  •  h) 

Vs*  1:3  '  53)* 

(124) 


•  ix)  Jp  *  Cs  •  iO 

•  iz)  ip  ^  (ip  •  iz)i,. 

•  iy)  ip  ^  (ip  •  iy)  i,- 

•!o)ip^ix^(i,-ix)ip-o^(V 

•  b\i  X  i  +(l  .  l\l  X  B  +/l  • 

•  B  \l  X  1„  +fl  .  i\l  X  B„  +^1  • 

~oj~p  yq  ^yj~p  ^0  ^~p 


1  \l  xB  +/l  •  b\l  xl 
~xJ-.q  -0  Up  “cjtq  tx 

5z)iq='?o+(lp’ 

iy)V!o*(ip-  !^q»iy 


Part  of  the  expression  (124)  for  the  magnetic  field  gradi¬ 
ent  depends  upon  the  vectors  Bx  which  are  related  to  the  gradl- 
P 

ent  Qpq  of  the  magnetic  field  in  the  absence  of  the  sphere.  It 
is  useful  to  obtain  Qpq  both  for  the  sake  of  comparison  and  for 
the  purpose  of  evaluating  (124)  in  more  detail. 
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As  in  Section  IV,  quantities  in  K  (transform)  space  will  be 
labeled  with  a  caret.  Thus,  the  vectors  6^  are  defined  as 
Fourier  components  of  derivatives  of  the  magnetic  field  even 
though  the  6^^  are  defined  only  at  the  center  of  the  sphere. 

Employing  (78)  and  changing  the  coordinates  from  C,  n,  C 
to  X,  y,  z  we  have 


B  -5^  6 

®Pz  ®Px’ 

Then,  from  (125)  and  the  fact  that 

V  •  Bp  =  0 


(125) 


it  follows  that 


and 


Now,  let 


Then 


-  i 


B 


Px 


K,  SBp 


^  -  i  K,  9Btt 

B  =  _  2  _ -JL 

3y  • 


B  = 


"  ^  ®Py  h 

3y  ^  y 


A  A  /s  ^  3  B 

5i-  S’  h  *  •'z  !>  ?3  ■  ^ 


From  (78c)  it  can  be  seen  that 


B, 


ay 


Py  W 


-  f  .[-1  K  .  B  +B  ^1 
-/_oo  ~  -0  oy  3yJ 


z-K|y-y' 


(126) 


(127) 


(128) 


^  (129) 
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It  follows  from  (127)*  (128),  and  (129)  that 


•D  op 

»  1  •  7(1  •  B  )  = _ 2. 

^pq  ip  ~q  ~P^  2K 


_£Yvit, 

^  ""  (130) 


where 


-V  'f  •  ip>  '5  •  5o>> 

<K  ■  B„)(K  •  1^)(K  •  1„) 

’’2  ■  -ipy  ",y  '!?  •  5o)  "  ~ -  V'5  '  ip>%’ 


’^3  “  ^py^qy®oy 


■^[ 


lpy'5  •  lq><5  •  ?o>  *  '5  •  lp>«  •  iq>®oy  > 


Ipy^oyt;  •  iq>. 


The  derivative  operations  In  (130)  can  be  carried  out  with 
the  aid  of  the  following  Identities: 


=  K[n(y’-y)  -  n(y“y*)] 
e  -  ‘'iy-y'l  -  -2K6(y-y')  + 


(131) 


_L^  e  "  Kly-y’l  ,  _2K6’(y-y')  -  K^e  "^1^"^’! 

•  Cn(y-y' )  -  n(y'-y)]» 

where  n(x)  Is  the  function  defined  by 


ri(x) 


y<  0 
y>  0 


n 


(132) 


-2KT 


n  n  " 


For  the  magnetic  field  gradient  Inside  the  sphere  (122) 
can  be  used.  The  result  will  depend  upon  the  moment  vectors 
given  by  (97) • 

Using  the  coordinate  system  with  the  origin  at  the  water 
surface  rather  than  at  the  center  of  the  sphere  then  (97)  can 
be  written 


0  00  00 


ii2 


^  f  f  f  p 


dx'dz'dy ' 


0  00  00 


uz 'dx' dz'dy * 


^+z  '^+(y'-y)^]^''^ 
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i^3 


'to(y '-y)dx*  dz*dy* 
Loo  J_oo  [x'  ‘^+z '  ‘^+(y  *-y 


(133) 


where  it  is  assumed  that  the  center  of  the  sphere  Is  at  the 
point  with  horizontal  coordinates  x  =  0,  z  =  0,  and  the  vertical 
coordinate  y  arbitrary.  By  considering  the  two  dimensional 
Fourier  transform  with  respect  to  x  and  z  it  can  be  shown  that 


Hi 


il2 


-K|y‘-y|  P 

e  w(K,y*  )dy'd‘^K, 


-K|y'-y|  P 

e  w(K,y*)dy'd‘^K, 


^3 


-K|y'-y I 

e  c-ri(y’-y)  +  n(y-y')]  u)(K,y  ’  )dy 'd^K 

mu  mu 


1  ±.  f"  f° 


/w 

(A) 


(K,y»)dy'd^K. 

^  mu 


(134) 


From  (124),  (127),  (128),  (129),  (134),  and  (8l)  it  is  found 
that 


pq 


■/> 


a  .3 

1  2  3y  3  jy'i  3yd 


“Kly '-y 1 

n 

(135) 


L  J  n  n2 


Where 

ay 

I'l  -  -5^  <5  •  ?oHVjy+K^A2,), 


49 


^3  - 


X  lay  r.  v 

■  5)«z<i2  •  ?>3<5  •  ?o>  -  [(Vly«z''2y)® 

*<S  •  ?o)'‘3y]’ 


cry 


2K 


^  (A3  •  K)(B^-  K)  +  [K^(A^- 


oy 


■2^  ^3yV* 


lay 


= 


25K 


[:^Kz^lx-Kx^lz)V^V2x-V2z>K,]. 


^2  =  i2^^-«z^3x-^Kx^3z) 


(136) 


The  p,q  dependence  of  G  is  determined  by  the  p,q  dependence 
of  the  A^  and  which  are  defined  in  connection  with  (124). 
Thus,  the  and  the  depend  upon  p  and  q,  although,  to  save 
space,  these  indices  are  not  exhibited  explicitly  here.  With 
the  aid  of  (13I)  G  can  also  be  written 


^  (L^  +  K2L3)e-''l>'’-yl2  4.„(y’)n„dy' 

n 

j^^e-Kly'-y|N2(y,)^  ^^4“ 

n  n 


i(L2  + 

n 

I  M,e~'^l^'~ylN^(y')^  ^nCyMUndy' 


-L3  2] 

n 


2^n(y)U  . 
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(137) 


VII. 


SPECTRA  OF  MAGNETIC  FIELD  GRADIENTS 


In  order  to  obtain  correlation  functions  ^^rs^*  which 

are  used  to  determine  the  spectral  characteristics  of  the  mag¬ 
netic  field  gradient,  it  is  convenient  to  write  (135)  in  the 
form 


[A. 


pq 


(y,y*  )y^ 

n 


d)  (y '  )U  +B 
^n^^  '  n  pq 


(y»y' ) ^  ^^Cy' 


)  U^]dy*, 


n 


0 


n 


(138) 


where 

Ap^(y,y') 


Ju 

2K 


t 


+  L. 


e-K|y‘-y| 


(139) 


and 

Bpq(y,y’)-  [M^+M2  e"^|y'"ylN^(y' ). 

In  (139)  the  quantities  and  depend  upon  the  indices  p,q. 

With  the  aid  of  Identities  derived  in  Appendix  E  of  Ref. 
1,  it  can  be  seen  that  the  correlation  averages  of  magnetic 
field  gradient  components  displaced  in  time  are  given  by 


-00  -00  (_  J  n 


{■ 


+|Ap^(y,y')B»g 


(y.y")^A5^<y.y")Bp,(y.y')]y; 

J  n 


kqCy.y  )Bf,<y.y")]2]  yCyMyty^ 

”  (140) 


J 
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where 


and 


V  =  i 
n  2 


in  (K)t  -in  (K)t 

,|;^(K)e  ""  +  e 


T  =  t*-t. 

If  the  Milder  hypothesis  [1]  is  used  in  (140),  this  becomes 

.  -/  r 

^ _ flo  _ fla 


“jTjL  V(y*y’)Aj3(y,y")^n;;(K)w^(K,T)4.n(y’)4>n(y 


")dy'dy" 


n 


A_(y,y’)B*g(y,y’')+A«^(y,y")B^g(y.y 


pq 


rs 


pq 


i 


23  (K)W^(K,T)(j>^(y’)<|)n(y”)dy'dy" 


n 


n  (141) 


where 


W„(K,t) 


n  •«. 


2TrK^  |_ 


r  in„(K)T  -in_(K)Tl 

I(K)e  +  I(-K)e  ""  J. 


(142) 


Now,  define 

Y^(y’,y”)  “23^n 


ly')  •J  Ajg(y,y”)  Y^(y’ >y'')dy". 


and 


B^g^(y,y’)  •  B»g(y,y”)  Y^(y ’ »y" )dy'’ . 


(143) 
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Then,  according  to  (1^1) 


=  yj| ■*■  4s^y»y’^j 

+Bpq(y.y')[A^^|^(y,y»)  +  B^°^y,y')]|  dy'. 


"rs 


and 


of  (143),  (131), 

fs? .  =1“)/. 

(1" 

^3rs^v 

(y',y)  -  LlJi 

"irs 
2K  . 

f  N2(y")e' 
■'-00 

^2rs 

2  . 

.-K|y”-y|  (y,^y„)jjy„ 


(1^5) 


'vy 


In  (145)  the  subscript  y  on  the  last  term  of  the  expression  for 

Apg^(y,y*)  Indicates,  as  usual,  differentiation  with  respect  to 

to  the  second  argument.  By  similar  means  the  relation  (l44) 
can  be  written 


-)/! 


',y') 


B^g^y,y')|dy' 


“  ^  y  e"^*^’"^*N2(y*)^A^J^(y,y»)+B^g^y,y’|ly' 


■(- 


.(V .  i>: 


.  ^ /VV‘'iy'-yiN2(y) , 

y 


A^^^y,y’) 


+B^s^y,y’)  dy' 


A^3^(y,y' (y,y' )|dy» 


-L 


4) 


3pq  rs 


(y,y)+B 


(2 

rs 


’'y*y>J-L||pq[ 


^(■1)  (y  =  > 

(146) 


,y')]c 

I  (y.y)l 


Again,  the  subscript  y’  In  the  last  term  of  (146)  indicates 
differentiation  with  respect  to  the  second  argument. 

By  comparing  (130)  with  (135)  and  using  (145)  and  (146)  an 
expression  like  (l46)  can  be  found  for  <Gpq(t)  G^3(t*)>.  The 
result  is,  in  fact,  (145)  and  (146)  with  and  Mg  set  equal  to 
zero,  the  quantities  and  replaced  by  ctvIqT^  and 

and  the  quantities  Lg  and  replaced  by  iay^Tg  and  lay^T^.  In 
this  replacement  the  complex  conjugates  in  (145),  indicated  by 
asterisks,  are  obviously  obtained  by  simply  changing  sign,  since 
the  are  all  real. 

Since  the  statistical  process  is  assumed  to  be  stationary 

actually  a  function  of  t  =  t’-t  in  its  time 
pq  rs  ys  ^ 

dependence.  If  <G„„(t)G^^(t * )>  is  integrated  over  the  two  dlmen- 

pq  rs 

sional  K  space  and  the  temporal  Fourier  transform,  with  respect 

to  T,  is  taken,  the  result  is  the  temporal  spectral  function 

4  (o) ,y  )  . 

pq;rs' 

If  this  process  is  applied  to  (146)  in  two  steps  then  the 
first  step,  integrating  over  K  space,  leads  to  a  sum  of  terms 
I^(T,y)  of  the  form 


I^(T,y)  -  /dy»  f  /"d^g /*  dy”  ^^^*7) 
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The  second  step,  taking  the  Fourier  transform  with  respect  to 
T,  proceeds  as  In  Appendix  E  of  Ref.  1. 

The  Integration  over  K  space  Is  done  In  polar  coordinates 
after  the  temporal  Fourier  transform  Is  taken.  That  Is,  the 
Integral  Is  taken  over  the  magnitude  K  and  the  direction  angle 
w  of  K.  The  Integral  over  K  Is  carried  out  explicitly  and  the 
angular  Integral  Is  over  the  Interval  (0,  -j)  after  the  change 
of  variable 


w  ®  a  +  B  -  IT 

from  w  to  6. 

When  this  procedure  Is  applied  to  (1^7)  and  the  derivation 
of  equation  (259)  In  Ref.  1  Is  Imitated,  an  expression  for  the 
magnetic  field  gradient  spectrum  analogous  to  (259)  in  Ref.  1 
will  result.  In  carrying  out  the  necessary  steps  the  large 
tow  speed  approximation  Introduced  In  Ref.  1  Is  essential. 

First,  It  Is  useful  to  consider  the  Integrals  J^(w,K,y), 
defined  by 

J^(w,K,y)  =^y'°H^(K,y,y',y")^  n^(t^(y’)<J)^(y")dy‘dy”,  (l48) 

-La, 

In  more  detail. 

It  Is  clear  from  (146)  that  the  possible  values  of  v  In 
(148)  are  0,  2  and  4.  Identities  (F-9),  (F-5),  and  (F-6)  of 
Appendix  F  In  Ref.  1  can  be  used  to  evaluate  the  series 

s^(y',y")  “2!^ 

n 

for  these  values  of  v. 

For  the  case  v  ■  0  the  expression  Is  simplified  somewhat 
because 


SgCy’.y")  ♦„(y')t„(y") 

n 


6(y»-y”) 
N2(y«)  ‘ 
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Thus, 


Jo(w, 


K.y)  -/*°  -  H^(K.y,y*  ,y')dy'. 


(150) 


For  the  other  values  of  v  the  calculation  Is  more  compli¬ 
cated;  (148)  becomes  for  v  ®  2,4: 

/*0  /-O 

J^(w.K,y)  H^(K,y,y',y'’)S^(y',y")dy»dy”,  (151) 

where,  according  to  (P-5)»  (F-6),  and  (F-8)  of  Ref.  1, 


S2(y' jy")  *  -  K  g(y* ,y”) 


Sn(y',y")  *  K^/*  N^(x)g(x,y'’)g(y' ,x)dx. 


(152) 


for  g(u,v)  defined  by 


g(u,v) 


e^^'^sinh  Ku  > 


(153) 


Inspection  of  (143)  and  (144)  shows  that 

H^(K,y,y' ,y')  =  Bpg(y,y’ )B»g(y,y ’ ) , 


(154) 


H2(K,y,y',y")  =  Apq(y,y*)B»g(y,y'’)+Bpq(y,y’)A»g(y,y''), 
Hi,(K,y,y' ,y”)  -  A  (y,y’)A»g(y,y"). 


,-Kly-y' 


According  to  (131)  and  (139), 

A(y,y')  -  -^1  (i^+A3)  +  K(L2+K^^)Cn(y’-y)-n(y-y’)]| 

-L3  6(y-y»)  -L^  «'(y-y») 


anu 

B(y,y')  -  N^(y')  I  M^+  KM2[n(y’-y)-n(y-y’)]|e,'^|y"y’ 


(155) 
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where  subscripts  p,  q  and  r,  s  are  not  shown  but  are  understood 
to  be  determined  by  the  constants  and  M^. 


It  follows  from  (150),  (151),  and  the  derivation  given  for 
the  relation  (259)  of  Ref.  1  that  the  temporal  spectrum  of  the 
magnetic  field  gradient  Is  given  by 

Jnn.r.,(a+®-'^»K)[I(K,a+B-iT)+I(K,a+6)]  (156) 

pq,x  a 


♦  (u.y) 

pq;rs 


“i  K 


dB, 


where 


J(w,K)  =  Jq(w,K)  +  J2(w,K)  +  J^(w,K). 


(157) 


The  magnetic  field  gradient  temporal  spectrum  in  the  absence 
of  the  bubble  can  be  obtained  from  (156)  by  setting  the 
equal  to  zero  and  replacing  by  au^  T^,  by  iay^T2,  by 
oVqT^*  and  by  iay^T^j  In  (155)*  The  result  is 


pqjrs 


(u,y) 


J^.,_(cx+6-Tr,K,y)  Cl(K,o+B-n)+I(K,a+B)  ] 


(158) 


+  J^^^(a-B-Tr,K,y)[I(K,oi-B-w)+I(K,o-B)] 


dB, 


where 


pq 


;rs  ./7 


°  (y,y’)A^*(y,y")Sij(y’  ,y'')dy’dy” 


•pq '"rs 


and  where 


A^(y,y’)  -  +  iK(T2+K^Tj,)[n(y’-y)-n(y-y’)]|e“*^'y"y’ 

-  ov^l  T26(y-y’)  +  IT^fi' (y-y’ )|  • 
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Each  of  the  terms  In  (156)  can  be  reduced,  so  that  the 

magnetic  field  gradient  can  be  expressed  as  an 

pq  y  rs 

Integral  over  g  of  a  quantity  that  Involves  Just  one  more 
quadrature.  For  this  purpose  It  is  convenient  to  define  func¬ 
tions  f^(K,y)  by  means  of  this  quadrature: 

/Q 

x'^e^*N^(x)dx.  (159) 

.7 

The  first  term  in  (157)  is  the  quantity  J^(w,K,y)  given  by 
(150).  The  relations  (154)  and  (155)  are  needed  in  evaluating 
that  term,  as  well  as  the  others.  It  is  found  that 


-T -  H  (K,y,y',y')dy‘ 

N^(y‘)  °  " 


N^(y*)e 


-2K|y-y' 


dy  • 


+  C 


o2 


^^ol'^^o2^®^^^/*°  e"^^y’N^(y')dy' 
"V  (A-2) 


+  < Co 1-^02 e^^y’N^(y')dy' 

%/  *00 

^^ol'^^o2^®^^^^o^"2K,y)  +  (C^^-C^2>®"^^^t^o^2^*""^"^o^2K,y)], 


where 


and 


-01  ■  i  "lpq”frs  +  ”2pq"?rs> 


=02  *  W  '"lpq«irs*''2pq"!rs>- 


(160) 


(161) 
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For  evaluating  the  second  term  in  (157)  It  Is  convenient 
to  define  a  f..ictlon  P(y',y)  by 


F(y 


.y)  =/ 


S2(y’,y")  A(y,y")dy”. 


(162) 


The  quantities  appearing  in  (l62)  are  given  by  (152),  (153), 
and  (155). 

Then,  because  of  (15^)»  it  is  found  that 
0  ^0  -  0 

jLXo  j  f®pq(y»y')^(y’»y)+Bre(y»y')Ppqfy'»y)]dy' 

(163) 

The  symmetry  of  the  function  S2(y*»y'*)  Is  used  in  deriving  (l62). 

Further  evaluation  of  the  second  term  in  (157)  proceeds 
from  (162)  by  means  of  the  equation 

-J  B(y,y')F(y*,y)dy’  *  I  (^  +  f^)e"^^^^^e^V(y’)F(y',y)dy’ 


|(Y  -  e‘^V(y’)F(y',y)dy'.  (164) 


The  subscripts  p,  q,  r,  s  and  the  complex  conjugates  in  (l63) 
can  be  added  to  the  general  relation  (l64)  by  using  the  appro¬ 
priate  subscripts  and  complex  conjugates  for  the  quantities  L^, 
^2*  l>y  Lij,  M^,  and  M2,  which  are  all  independent  of  the  inte¬ 
gration  variables  in  (162)  and  (163). 

The  evaluation  is  continued  by  proceeding  from  (162). 

Thus,  using  notation  the  meaning  of  which  should  be  obvious. 


p(yt^)  .  -kKL^  slnh  Ky  +  L|j  cosh  Ky)e^' +(021-022)6"'^ slnhKy 


.-Ky, 


slnh  Ky”  dy” 
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+(C2i+C22)e^^^'‘^^^/'V^'’slnh  Ky”  cSy’*] 

•'y 

=  -K  j-CL^  sinh  Ky  +  KL^  cosh  iQr)e^*+  ^  (C2^-C22)e"^e^^’  sinh  Ky' 

+  I  (C2i-C22)[  ^  +  y-yDe^^^’-y^  +  |(C2i+C22)e^^y^' ^ 

X  C-y  +  ^  j 

=  PQ(y)e^'  +  Pi(y)y’e^’,  (I65) 

where 

Q 

pQ(y)  =  KCCL^  -  )  sinh  Ky  +  KLj,  cosh  Ky  +  |  (C23^+C22)ye^ 

+  1  (C2i-<322)y®’^» 

PjCy)  =  - 

Cji  ■  S  (167) 

and 

C22  *  I  (^2  +  K^Lij) 

Similarly, 


P(y'  >  y)  •  -K  j-(L2+KL^^)e^slnh  Ky»  +  sinh  Ky»  [(023^-022)6"^^ ^e^*^"dy” 

«  dy«  j 


+  (C2i+022)e^/’^  dy"]  +  f' 

*^V  •«* 


e"^"slnh  Ky’ 


-  -K  j  -(L3+!a,^)e‘^  *y’  C  ^  (Cai-Caz 

♦  f(02i+Cjj)e‘«y'*y>[-y  +  ^  (l-e-2^')3  j 

-  P2(y)e*^'-P2(y)e"^  +P3(y)y’e"^’, 


)+(023^+022)(y’-y)] 

(168) 
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where 


F2(y)  =  I  [L^  +  KL^  -  I  +  C^^)y2e^ 


(169) 


F3(y)  *  I  +  C22)e^^. 

Inserting  (I65)  and  (168)  into  (l6^) 
y*  B(y,y')F(y',y)dy'  =  |  (^  "^  M2)e"^|  F^(y)[f^(2K,-o)-f^(2K,y)] 

+  F3^(y)[f3^(2K,-»)-f^(2K,y)]|  (I70) 

I  F2(y)[f^(o,y)-f^(-2K.y)]+P3(y)f^(-2K,y)j. 

The  desired  expression  (I63)  is  obtained  from  two  applications 
of  (170). 

The  subscripts  on  B  In  (163)>  as  well  as  the  complex  con¬ 
jugation,  are  obtained  by  applying  the  same  subscripts,  and 
complex  conjugation,  to  the  quantities  and  M2  appearing 
explicitly  In  (170).  The  same  rule  applies  to  A  with  subscripts 
and  complex  conjugation  applied  to  L^,  L^,  C2^,  and  0^2  appear¬ 
ing  explicitly  In  (I66)  and  (I69). 

In  the  last  term  of  (157)  the  expression  to  be  evaluated 
Is,  because  of  (151),  (152),  and  (154), 

rr  H4(K,y,y',y")Si,(y’,y")dy”  -  l^{x)G^liyyX)G!^{y,x)6x,  (17I) 


where 


r,x)  »/  A(y, 
.00 


y’ )g(y’.x)dy’ . 


(172) 


In  (171)  the  subscripts  on  Q(y,x)  and  the  complex  conjugation 
actually  occur  through  A(y,y')  In  (172). 
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With  the  aid  of  (155)  the  expression  (172)  can  be  written 


(173) 


G(y,x)  »  -L2g(y,x)  -  L2jgy(y,x) 

+  (C2i-C22)e"^y'^e^^’g(y’,x)dy*  +  (C2i+C22)e^y^°e~^y’g(y',x)dy’ 


From  (153)  it  follows  that 
L, 


G(y  >  x)  =  -  (-j^  sinh  Ky  +  ^  cosh  Ky)e^  +  ^  (C2^-C22)e“^[sinh  Kx 

(174) 


xA^Ky* 
•^x 

+  i<A^  +  Aj)  f  e'^’sinh  Ky'dy* 


+e*“‘(  sinh  Ky*  dy*] 

X 


*  GQ(y)e^  +  G^(y)xe^, 


where 


G„(y)  -  i  -  ye'^)  +  i  (c„-c„)(2M 


2K  ''"21  '"22"'  K 


-ye-*^) 


,  sinh  Ky 

-L3  ^  -  L 


and 


rjij  cosh  IQr 

Gi(y)  =  ^  (C23^-C22)e"^. 


(175) 


It  also  follows  that 
L 


G(y<x)  ■  -(-j^  +  Lj^)e^  sinh  Kx  +  ^  ^^2l”^22^®”^  ^ 


+  1 


[sinh 


e“^'sihh  Ky’dy’]  (176) 


G2(y)e'^  -  G2(y)e-^  +  G3(y)xe’^, 
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(177) 


where 


Q,(y)  =4  ^  -  La  -  F  (0„+C„)y]e’^ 


K  "'21  "22' 


and 


G-Cy)  =  -  (C2i+C22)e^. 


It  follows  from  (171),  (174),  and  (176)  that 


/T 

oo^'.oo 


Hi,(K,y,y',y")Sij(y*,y”)dy’dy" 


=  k'*//n2(x)  [G^j^(y)eK^+G^pq(y)xe'^][G2^(y)e^+Gj^(y)xe'^]dx 

»  K^l  H^(y)Cfo(2K,-«)-f^(2K,y)]+H^(y)[f^(2K,-«.)-f^(2K,y)] 
+H2(y)[f2(2K,->)-f2(2K,y)]+H^^(y)f^(-2K,y)+H^(y)f^(o,y) 
■^”ol^y^^o^^^»y^'^%o^y^^l^°*y^‘^”ll^y^^l^2K,y)+H2^(y)f2(2K,y)|  ,  (178) 

where 

«o<y>'°opq<y>°;rs<5'>’«l<y>  =  °lpq<y>°Srs'y>*°oixj'y>°irs<y>> 
•*2<y>*°lpq<i'>“!r5<5'>>"o-l<5'>’°2pq<y>°lre<y>>“oo'y>‘-“2pq'^'°§rs<y>* 
«ol'J'>-°2pq<y>“5rs<y>>«lo<y>  ■  -°2p,<y>°!rs<y>-°3pq'yW|rs'y>> 
“ll'y>'02pq<y>°|rs'y>«3pq'y>°lrs'y>’“2l'='>  "  °3M<y>°!rs<y>-  <179) 


When  there  is  no  sphere,  the  magnetic  field  gradient  is 
that  due  to  Internal  waves  for  the  case  in  Ref.  1  where  the 
field  is  observed  below  the  sea  surface.  When  the  gradient 
contains  a  vertical  (y)  derivative  it  will  have  a  discontinuity 
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w 


at  the  sea  surface,  but  If  the  derivative  is  in  the  horizontal 
(x  or  z)  direction  it  should  be  continuous  at  the  surface. 

Thus,  the  quantity  given  here  for  the  case  of  no  sphere 

should  agree  with  (259)  in  Ref.  1  when  both  are  evaluated  at 
y  =  0  if  the  assumptions  that  were  made  in  Ref.  1  in  deriving 
(259)  are  also  made  here.  These  assumptions  include  the  Milder 
hypothesis  and  a  large  tow  speed,  both  of  which  have  already 
been  adopted  here.  In  addition  the  assumption  of  an  isotropic 
internal  wave  excitation  function  with  power  law  dependence  on 
wave  number  was  also  used  in  Ref.  1.  That  is,  it  was  assumed 
that 


I(K)  =  CK"P 


For  the  case  of  no  sphere,  as  already  observed  in  con¬ 
nection  with  (157),  J  reduces  to  J|^,  and  at  y  =  0,  according 
to  (175)>  reduces  to 


2 


V^^21pq"^22pq^^^21rs"^22rs^^2^^^»"”^ 


K^c 


I- 


K 


KT, 


3pq 


+  ^3rs*iT2rs 


) 


f2(2K,-«). 

(180) 


For  the  last  expression  in  (l80)  the  rules  for  substituting 
T^  in  the  case  of  no  bubble  were  used  in  (167). 

According  to  the  definition  of  the  Tjj^  in  connection  with 
(130)  and  the  dependence,  given  by  (113)  in  Ref.  1,  of  on 
<|)p,  and  by  (120)  in  Ref.  1  for  the  Ip  vectors  on  a,  it  is  a 
straightforward  matter  to  verify  that 

*^113  “  *^213  “  B^Kcos(w-a)sln(w-a)cos  w  cos  (Ifil) 

“  "®o  (w-a)  sin  (w-a)  sin  i|)p. 
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If  (181)  is  substituted  Into  (I80)  it  is  found  that 


2a^y% 

‘T)-' 


where  defined  by  equation  (l82f)  in  Ref.  1 

Finally,  if  (182)  is  used  in  (I58)  along  with  (179) 
will  be  found  that  the  result  agrees  with  (259)  of  Ref.  1 
order  to  complete  the  verification  it  is  necessary  to  use 
to  identify  f2(2K,-«). 


(182) 


it 

.  In 
(159) 
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PROOF  OF  THE  DIRICHLET  GREEN'S  FUNCTION  APPROXIMATION 


The  notation  to  be  used  in  the  following  is  consistent 
with  that  of  the  preceding  text.  Namely,  Gg(r,r*)  is  the  exte¬ 
rior  Green's  function  for  S  (G„  vanishes  on  S)  and  G^  (r,r' )  is 
the  Green's  function  for  Vp  (Gp  vanishes  on  P).  In  fact,  with 
GQ(r,r')  the  free  space  Green's  function  given  by 


G 


o 


9 


G^CrjT')  -  G^(r,r') 


and 


4Wa^-2a^rr  *  cosy+r^r ' 

In  addition,  operators  Dg  and  Dp  are  defined  by 
D.G.-/*  G.(r",r')ds" 

S  i  y  — IF’ -  ^  ^ 


(A-1) 


(A-2) 


and 


D.,G^»  /"  G^  (r",r' )d5".  (A-3) 

^  ^  J  —dP' -  ^  ''  ' 

S 

According  to  Green's  theorem  (A-2)  implies  that  DgG^^  is  equal 
to  G^  on  S  and  satisfies  Laplace's  equation  outside  of  S. 
According  to  (A-3)  DpG^  is  equal  to  DgG^^  on  P  and  satisfies 
Laplace's  equation  in  Vpg. 


A-3 


nacnoao  njNMiot  iumd 


Let 


“I'r-r')  =  ODs'r-r'>  -  Ss'r-r')-  (a-*!) 

Then  vanishes  on  P  but  not  on  S.  However,  ^^^^ishes 

on  S  but  not  on  P.  Further,  G^-D^G^  +  DpG^  vanishes  on  P  but 
not  on  S,  and  +  DpG^-D^DpG^  vanishes  on  S  but  not  on  P. 

Continuing  the  Iteration  ad  infinitum  leads  to  a  Green's  func¬ 
tion  Gp(r,r*)  defined  by 

00 

Gjj(r,r’)  =  (I-D3)  ^  Dp*^  =  (I-D3)  (l-Dp)"^^  ^  (A-5) 

n=o 

if  the  series  converges.  The  limiting  function  vanishes  on 
both  S  and  P,  as  required  for  the  Dlrichlet  Green's  function 
appropriate  to  the  region  Vpg  bounded  by  P  and  S. 

Convergence  of  the  series  in  (A-5)  occurs  in  the  sense  of 
any  norm  for  which 

||Dp|l.<l.  (A-6) 

For  any  function  f(r") 


A-i» 


so  that 


(Pi-7) 


if  the  norm  of  a  function  is  the  maximum  value  of  the  function 
on  S.  Convergence  is  then  guaranteed  as  long  as 


a 

D 


< 


1 

3* 


(A-8) 


In  forming  the  series  that  converges  to  G(r,r')  each  step 
adds  a  term  one  order  higher  in  In  fact,  the  series  begins 
with  G^,  the  free  space  Green's  function  and  G^  ■  (l-Dg)G^. 
Thus,  the  lowest  order  approximation  to  Gp(r,r*)  Is  G^(r,r') 
which  is  also  given  by  (A-4). 

By  definition. 


AG  *  G»>**•G^T^ 
DP. 


Hence, 


so  that, 
by  (60), 


AG  ~ 

according  to  (A-1)  and  (A-4), 
as  asserted  earlier. 


AG  Is  given  approximately 
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In  the  following  it  will  be  shown  that  the  partial  electric 
field  contribution  J^.  to  the  current  given  by  (4?)  results  in  a 
contribution  of  higher  order  than  a  to  the  magnetic  field  in¬ 
side  the  sphere.  Then  it  will  follow,  once  again,  that  only 
the  term  in  (47)  need  be  considered. 

The  potential  4>’  whose  gradient  determines  Jg  satisfies 
Laplace's  equation  inside  Vg  and  the  boundary  condition 


ii'l 

3r  'r=a 


ir 


r»a 


(B-1) 


on  S.  From  (B-1),  (53),  and  (64)  it  follows  that,  to  terms  of 
order  less  than  a. 


(i.  •  M)  <1_  •  B„) 


r=a 


a 


+  i  •  P, 


(B-2) 


where 


?  -  <5o  • 


a 


(B-3) 


In  obtaining  (B-2)  terms  of  order  ^  have  been  dropped  as  usual. 
Now,  for  any  two  constant  vectors  and  kg  the  quantity 
given  by 

-  'Ir  • 

is  a  surface  spherical  harmonic  of  the  second  degree  [33*  Also, 
for  any  constant  vector  k  the  quantity  i;  given  by 

V  -  .  k  (B-5) 

1 

B-3 


Is  a  surface  spherical  harmonic  of  degree  one,  and,  of  course, 
any  constant  ^  ^  is  a  surface  spherical  harmonic  of  degree  zero, 

It  follows  from  (B-2),  (B-4),  and  (B-5),  for 


Zo((|),0)  * 


0 

li  * 


Zi((|.,0)  »  Ip  *  P, 


(B-6) 


and 


ZjC^.e)  -  |[(lr  •  vXir  •  ?o)-  Xh  •  ?o)]- 


that 


.n+2 


n+1  ^a^****®^ 


(B-7) 


In  (B-7)  the  quantities  (ij),0)  are  surface  spherical  harmonics 
of  degree  n. 

Prom  (B-7)  it  follows  that 
2 


7<^' 


so  that 


ip  X  V*' 


Thus, 


[i. 


7oZ„(^,e)  ■ 

IJ 

n+1 

(i 

if"  ir  ='  . 

V 

Fj  n+l 

n-1  ' 

/ 

,  X  7<fr» 

J  r«a  tl  i2* 

(B-8) 


B-M 


(B-9) 


where 


and 


X  V 


n 


P) 


X  V 


a 


la 


Prom 


V(r  •?)»?=  V(ri  •  P)*  (1^  *  p)i  +rV(l  •  p) 
■  <Jr-  ?>ir  *  ''a  <ir  '  V 

It  follows  that 

ir  *  ^ir  *  ?>  “  ir  *  ?• 

Since  (B-10)  is  an  Identity  ths.c  holds  for  an  arbitrary  con¬ 
stant  vector  P  It  can  be  applied  to  I,,  with  the  result  that 


Ja  ■  -  35  »  C'lr  ■  ?o’''ll<ir  ’  ">+<ir  '  “'’n'ir  ' 


■  -  3a  Jr  *  C'lr  '  ?<,>“  *  'Jr  ' 

It  Is  useful  to  define 
-  "  ?o  *  i!»  ■  V  X  u,  X,  «  V  X  . 


Then 


X-  X  V 

P  -  — g  ,  B 


X-  X  V 

C  a« 


-0 


V  -  V 

Prom  (B-11)  and  (B-12)  It  follows  that 

X 


(B-11) 


(B-12) 


(B-13) 


h--ii  [(ir  ■  ?o)(lr  •  {ir  '  ?o)(‘r  '  h)^ 

*  <ir  •  HHlr  '  X’  7  '  (Ip  '  v)(l,  ■  Xj)  .  (B-U) 


From  the  form  (B-4)  It  can  be  seen  that  each  cartesian 
component  of  I2  Is  a  sum  of  a  spherical  harmonic  of  degree  two 
and  a  spherical  harmonic  of  degree  zero  (constant).  A  simple 
calculation  shows  that  the  spherical  harmonic  of  degree  zero, 
in  fact,  vanishes.  Thus,  I.  is  a  vector  whose  cartesian  com- 
ponents  are  each  spherical  harmonics  of  degree  two.  Because 
of  (B-10)  it  is  apparent  that  I.  is  a  vector  whose  cartesian 

J. 

components  are  each  spherical  harmonics  of  degree  one. 

Then,  as  indicated  by  (52),  (51)  becomes 


?E'r=a  ■  lir 


=  au^a 


r/'n 

n*l 

j■I^((l),e)  ^  l2($,0)j^ 


(cos  y)  i^((|)*,e*)dn 


(B-15) 


by  virtue  of  an  Identity  given  in  Ref.  3»  p.  137.  It  follows 
from  the  definition  of  in  (B-9),  from  (B-10),  and  from  (B-11) 
that 


Je  U  ir  ^  ^  I5  ir  * 


[<! 


•  B  )m  +  (1 

I*  0 '  ^ 


(B-16) 


In  accordance  with  (38)  and  the  earlier  remark,  in  connec¬ 
tion  with  (48),  concerning  the  meaning  of  I^,  the  only  contri¬ 
bution  of  Jg  to  the  magnetic  field  in  Vg  is  proportional  to  the 
radial  component  of  I-,.  It  follows  from  (B-I6)  that  this  con- 
tributlon  is  zero  up  to  terms  of  order  less  than  a  ;  thus,  as 
promised,  the  contribution  of  J_.  to  the  magnetic  field  can  be 
Ignored  in  calculating  the  field  to  the  first  order  in  a. 
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Formulas  needed  for  calculating  the  magnetic  field  gradi¬ 
ent  spectrum  appear  In  the  main  text.  In  fact.  It  Is  only 
necessary  to  apply  the  following  equations,  considered  In  the 
following  order  rather  than  the  order  In  which  they  were  orig¬ 
inally  Introduced,  to  establish  the  necessary  computational 
logic:  (156),  (157),  (151),  (160),  (l6l),  (159),  (163),  (170), 
(166),  (169),  (167),  (178),  (179),  (175),  (177),  (135),  (12^4). 

The  key  parameters  In  the  calculation  are  given  by  equa¬ 
tions  (136),  reproduced  here  for  convenience: 


L.,  »  - 


2K 


ay 


2K 

lay 


1  lay  r.  v 

•  ?o>  - 

«3  •  E)3Boy 


2K 


°  <*3  •  E>®oy- 


"  "3y®oy 

(136) 


ay 

”i  "  ^  C(K2r^^-K^r^^)K^+(K2r2x-Kx^2z^^z^* 


”2  “  ^^^•^z^3x^^x^3z^’ 
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A 


For  purposes  of  comparison  and  interpretation  the  various  quan¬ 
tities,  (136)  can  be  expressed  in  terms  of  geometrical  and  phys¬ 
ical  parameters  that  were  Introduced  in  [1].  Each  of  those 
quantities  is  a  function  of  two  Indices  that  are  determined  by 
the  subscripts  on  the  symbol  for  the  magnetic  field  gradient 
spectrum.  The  indices  of  Interest  are  the  pairs  12,  13,  23; 
for  these  index  pairs  the  following  quantities  appearing  in 
(136)  are  defined: 


V12  "  Vl3^  ^  sln(w-2a),  K)2i  =  0. 

(A2-  K)i2  =  0,  (A2-K)3^3»  K  cos(w-2a),  (A2-K)23“ 

(A^*  K)^2  **  K  cos(w-a),  0,  (A^*  K)23  ®  K  sin(w-a), 

^lyl2  “  ^Iyl3  *  ^ly23  “  ’ 

^2yl2  ""  sin  ^2yl3  ”  ^2y23  ”  cos  o(>  (C“l) 


^3yl2  “  ^3yi3 


^3y23  ~ 


1x12 


^1x13  ’  "®o  2a,  r^^23  “ 


^lzl2  "  -SB^cos^pCOsa,  -B^sin(|>pSin2a,  2B^cos(|)pSina, 

^2x12  “  B^cos<|ipCOSo,  =  -B^sin<|)jjSln2a,  r2j^23“  “B^cos^pSlna, 

^2zl2  “  -B^cos^pSlna,  B^sln(^jjCos2a,  r2223"  -B^cos^jjcoso, 

^3x12  “  “SB^sln^pSlna,  ■  -BQCOS(fjjCos2o,  r3j^23“  -2BQSln<^jjCosa, 

^3zl2  “  2B^sln(|>jjCOSo,  - -B^cos({»jjSin2a,  ^2z23  "  -SB^sin^DSin®* 


C-iJ 


where 


Equations  (C-i)  used  in  (136)  define  all  of  the  and 
corresponding  to  the  subscript  pairs  12,  13,  23.  These  are  the 
parameters  to  be  used  when  the  sphere  exists.  In  the  absence 
of  the  sphere  the  are  set  equal  to  zero,  and  are  re¬ 
placed  by  and  and  and  are  replaced  by  l(JyoT2 

and  lay^T|j,  where 

■^Ivl  “  '^lv3  “  ^12p  “  '^112  “  -K^B^cos(w-a)cos  w  cos  (|)p, 

2 

(w-a)  cos  w  cos  (j)^,  ^221  “  "^223  “ 

2 

T211  *  KB^cos  (w-a)cos  w  cos  <|>p,  ^2^2  “  -KB^cos  (w-a)sln  4>p, 

T213  =  2(w-a)cos  w  cos  T222  “  -KB^cos  w  cos  ((tjj, 

T231  ■  2(w-a)cos  w  cos  (^p,  T222  “  -KB^sln  (w-a)sln  (j>p, 

2 

T233  *  KB^  sin  (w-a)cos  w  cos  ({)p>  T2j^2  “  *^332  “ 

2 

®  -Bq  cos  (w-a)sln  (|)p,  “  "®o  2(w-a)sln  (J>p, 

1^21  ®  "B^  cos(w-a)cos  w  cos  (J)p,  "^222  “  ®o  *^D* 

T323  ®  -B^  sln(w-a)cos  w  cos  4ip, 

1  2 
T221  *  -  2®o  ^  (w-a)  sin  (j)^,  “  "B^  sin  (w-a)sln  4)p.  (C-2) 


The  quantities  or,  in  the  case  of  no  sphere, 

the  are  used  to  define  the  quantities  and  Co2»  given 

by  (161),  the  quantities  C22  and  C22  given  by  (167),  the 
functions  PQ(y)  and  Fj^(y)  given  by  (I66),  the  functions  F2(y)  and 
F2(y)  given  by  (169),  the  functions  GQ(y)  and  G^(y)  given  by 
(175),  and  the  functions  G2(y)  and  G^^y)  given  by  (177).  For 
convenience  these  equations  are  all  repeated,  as  follows: 


'ol 


1  f  1 

t  '"T 


M.  M?  + 
Ipq  Irs 


2pq  2rs' 


=02  ■  W 


(161) 


(167) 


^21  “  2K  ■'■ 

C22  =  I  (^2 

n 

F^Cy)  =  KCCL^  -  -|i  )  sinh  Ky  +  KL^  cosh  Ky  +  |  (C2^+C22)ye^^ 

+  I  (C2i-C22)ye"^^,  (166) 

_  F^(y)  =  - 

F^iy)  =  I  CL3  +  -  I  ^23^  +  (C23_  +  C22)y]e^^  (169) 

F3(y)^  =  I  +  C22)e^^. 

G,(y)  =  i  (C2,k:22)(^^  -  ye^)  .  i  (C2,-C22)(^^  -  ye-^) 

-L^  ^  cosh  Ky 

Gi(y)  =  ^  (C23_-C22)e"^. 

GsCy)  =  ^  [^  -  ^  -  .  1  (C2^+C22)y]e^ 

^  ~  2K  (G2i+G22^®^’ 

In  addition.  It  is  convenient  to  define 

1 

^  -I  ^2^* 

For  the  case  in  which 

I(K)  -  CK"P 


(175) 


(177) 


(C-3) 
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equation  (156)  for  the  magnetic  field  gradient  spectrum  can  be 
written 


^  ^  /  2  ^-(P  +  1 

vy;rs  m 


dB 


where 


de 

(C-l|) 


^(°2vu°!rs  -  °ovp  °Jr3>fo<2>'.y) 

*  °ovu°Jr3  -  «=2vB°!i-s*03vu°lr3>fl'°*i'> 

^'°2vu°W°3vM«lrs  -  “ivv^Srs  '  aovy^Jrs'^l'^K.y 

*<«3vy°3rs  '  °lvy0!r3>f2'2>'>y>  *  “ivv^iry 


and 


■'iilrs'"’*'’  -[<=01  ^  =02>*"‘'’'  -  <I>Jy  ”^3  ^  V  "bw fo' 

{<=02  -  =ol>®'"^  -  <=vy  ''Srs  *  =J3  ’’ovy>«'^’']=o 
<=ol  -  =o2)«‘"'='  *  <=vy  ^Sra  "  “Js  ''ovy>'’^‘'o 


(D»  Pi  +  D  P,  )e 
'  vw  3rs  rs  3vy 


-Ky  f  (_2K, 

J 


y) 
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00  )j 

-2K,y) 

(2K,y) 

(2K,--) 


-['“vM  ''!rs  *  DJ3  F,„^)e-‘'yjf,(2K,y) 

"irs  *  “?=  (C-5) 

expanded  in  terms  of  the  moments  defined  by  (159),  l.e.,  the 
f unct ions 


f;^(K,y)  =/  e^*N^(x)dx. 

''y 

To  save  space.  In  (C-5)  the  y  dependence  of  the  functions 
Pj^(y)  and  Gj^(y)  Is  not  Indicated.  In  (C-4)  the  integrand  Is 
separated  Into  the  two  functions  and  because  in  the 

absence  of  the  sphere  vanishes  leaving  only  suitably 

modified  by  the  rule  requiring  the  replacement  of  the  quantities 
with  and  the  omission  of  the 
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